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SUMMARY 

This report describes a FORTRAN coded computer program and method to 
analyze the statistics of a transient 'capture* manuever for a single axis 
on-off controlled upper stage of a rocket powered vehicle. It uses a simple 
integration of the equations of motion to compute the time history of the 
transient maneuver. The Monte Carlo technique is employed to determine the 
statistics of maximum attitude error overshoot, time to maximum overshoot and 
total impulse expended. Booster induced disturbances, initial conditions, and 
aerodynamic disturbances are the main sources of excitement. These are 
treated as random variables. Certain vehicle and control system 
characteristics are also treated statistically. 

Control logic assumes a simple on-off attitude control system with a 
deadband. An optional limited integral of attitude error term is included in 
the control logic. 

Pitch, yaw and roll axis transient maneuvers can be evaluated by 
appropriate selection of input data. 

Running time for a Monte Carlo analysis of 199 samples is 5 seconds on 
a CDC Cyber 175 computer system. 
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1.0 INTRODUCTION 


Evaluation of reaction control system performance includes transient 
responses to events such as stage separation transients, rocket booster motor 
ignition, or handoff from a previously operating control system. Peak 
attitude error, time, and impulse experienced during these transients are a 
measure of the control system performance. A vehicle experiencing a rocket 
motor ignition while aerodynaiaic moments are still appreciable can diverge and 
lose control. Figure 1 presents a phase plane trajectory of a capture 
maneuver for the Scout launch vehicle. A successful and a divergent case is 
shown. The divergence results from aerodynamic disturbances. Probability of 
successfully capturing during these events is usually required in design 
analyses or for definition of permissable staging criteria. 

Designing to worst case criteria is sometimes used but can result in 
severe penalties in the system performance. A statistical approach is usually 
preferred to gain a sound basis for control system design with qualitative 
results. Non-linearities in the disturbances and control system behavior 
result in a difficult statistical prediction unless a Monte Carlo analysis is 
used. The Monte Carlo technique is employed in the ’capture’ transient 
computer analysis routine described herein. This method has been used for 15 
years on the NASA SCOUT launch vehicle upper stage control systems. 

The accuracy of this method is only as good as the validity of the 
assumptions as applied to the particular case and the accuracy of the input 
data statistics. It is suggested that sensitivities to variation in the 
random variables be evaluated for a given system. Occassionally a modified 
Murphy’s Law prevails, i.e. the most significant random variable is the one 
for which the least accurate information exists. 



2.0 METHODOLOGY 


Prediction of the staging capture maneuver statistics presented herein 
makes assumptions which may or may not apply to a particular control system. 
The method is for a single degree-of-freedom (pitch, yaw, or roll) time 
response of a launch vehicle upper stage using an on-off type reaction control 
system. Assumptions and equations for the control system, vehicle, and 
forcing functions are described in this Section. Statistical analysis 
methodology is also presented. 


2.1 Assumptions 

A block diagram of the general control system logic is presented in 
Figure 2. Some of the assumptions made in this method are: 

. single axis control unaffected by other axes motion or control, 

. mass properties do not vary during capture maneuver, 

. velocity and wind velocity can be defined by a linear variation 
with time (two point curve), 

. control motor thrust is constant level and it responds to 'on' 
and 'off' commands with separate time delays, 

. gyro and other guidance system component responses can be 
represented by a time delay, 

. aerodynamic normal force coefficient can be expressed as an odd 
function of angle of attack including linear and cubic terms, 

. aerod 3 rnamic center longitudinal location is an even function of 
angle of attack and can be expressed as a constant with linear 
variation with the absolute value of angle of attack, 

. vehicle bending effects are represented as a simple control motor 
switching delay time based on a random phase of one bending cycle, 

. initial attitude error includes (l) a bias value, (2) a function 
of initial angle of attack, and (3) a random sinusoidal component 
with an amplitude having a Log-Normal Distribution Function, 

. the vehicle is statically balanced and is aerodynamically 
symmetric. 
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2.2 Pitch Plane Equations of Motion 


2.2.1 Angular Accelerations 


Angular equations of motion written in terms of the accelerations are 
presented below. See Figure 3 for sign convention. The total angular 
acceleration includes control, aerodynamic and booster thrust misalignment 
disturbances. 


(2-1) 0= + 0’aero+^r Cdeg/sec2) 

«• 

where 6q is the control angular acceleration 

(2-2) ^c= 57.3 (Ft - Fb)[(.Xc - Xeglcos\^+ Zq sin / 12 ly (deg/sec2) 

^aero aerodynamic disturbing acceleration 

(2-3) ^aero = 57.3 CNSqCXgg - Xcp) /.12 ly (deg/sec2) 

•« 

d^i.3 the rocket booster motor disturbing acceleration 
(2-4) 0 ^ = (Xr - Xcgl / 12 ly (deg/sec2) 


2.2.2 Aerodynamic Coefficients 


Aerodynamic coefficients are assumed to be the following function of 


angle of 

attack 


(2-5) 

CrS = (CujjS)a=o“+ 

(ft2) 

(2-6) 

1 „ 1 

^cp - Xcpti=o 5a 1 1 

(inches) 


This provides an approximate form to the supersonic aerodynamic 
coefficients of a symmetric missile type body at angles of attack over a plus 
or minus 10 degree range, lypical aerodynamic coefficient variation with 
angle of attack are shown in Figure 4. If it is not available the non-linear 
coefficients can be approximated by the method of Reference 1. 
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2.2.3 Angle of Attack 


Angle of attack is computed from the trajecfcory flight path deviation, 
attitude error, wind, velocity, and path angle by the following relationship 

(2-7) o “ ^ ^wp (degrees) 

Pitch attitude error is obtained from the attitude minus the commanded pitch 
program, 


(2-8) 

de = d - dc^ 

(degrees) 

Wind and trajectoiy parameters are assumed to be linear functions 
following ignition, i.e. 

(2-9) 

y = y„ + Set 

nominal flight path angle (deg) 

(2-10) 

V = Vo + vt 

velocity (ft/sec) 

(2-11) 

= ^VPo t W 

wind velocity (ft /sec) 

(2-12) 

“vp ”®WPo t “wpt 

angle of attack due to wind (deg) 


Angle of attack due to a horizontal wind component in the pitch plane (V^p) 
is, 

^degrees) 


. (degrees) 


Conditions 


(2-13) = 57.3 


'vp 


sin Y 


V + Vyp cos y 


Angle of attack can be expressed as 
(2-14) a = 0 + - Ye * (“wp - 


2.2.4 Initial 


Initial conditions for a vehicle at an upper stage ignition will depend 
upon (l) the behavior of the previous stage and (2) the 'tip off' or change in 
characteristics due to an impulsive type ignition or separation disturbance. 
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If the previous stage has a proportional control system and aerodynamic 
disturbances are significant, as in the case of the Scout second stage, the 
initial conditions will be correlated with the trajectory deviations and wind 
velocity. It is assumed that the initial attitude error can be expressed as a 
combination of static trim and an oscillatory motion as shown in the sketch 



(2-15) ^trim+ ^amp sin T] (degrees) 

where tj is a random phase angle for the oscillatory motion 

Trim conditions of the prior stage include a nominal bias and a bias 
due to the steady state angle of attack due to wind and trajectory errors. 

(2-16) 0trim = ^e^ias (degrees) 

where 6*ebias previous stage trim attitude error due to vehicle 

asjrmmetries (degrees) 

Kfl is the sensitivity of previous stage attitude error to angle of 
attack 

Initial angular rate must include the total from previous stage motion 
and tipoff disturbances. 


2 . 2.5 Control Logic 

Control logic is according to the block diagram of Figure 2. The error 
signal is given by 

(2-17) 6p = Oq + Kr/Kd 0 + GjK^e<^^j (degrees) 

where Kr/Kr is a constant rate to displacement gain ratio. The integral 
gain term can be delayed and limited. Integration of the attitude error 
starts at ignition, or time zero. The integral gain (Kj) is a constant but 
can be set to zero for an initial time period (tgj) following ignition. 

This helps to alleviate the initial transient capture. Limiting of the total 
integral gain term can also be used. This usually is desirable to make the 
control system dynamically stable. The integral gain function can be 
expressed by the following equations 
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(2-18) 


t ^ 0 


For, 

(2-19) 

The integral gain term can be deleted completely by input option selection. 

Control motor firing logic based on the error signal generated by 
equation (2-17) and the deadband halfwidth (d) is as follows: 




Fb 

Frp 

(2-20) 

Sp > d 

ON, 

OFF 


-d < 6p < d 

OFF, 

OFF 


6p < -d 

OFF, 

ON 


Equation (2-20) defines the nominal switch points (also see Figure l) . 
Ifysteresis in the switching is treated by narrowing the switch off deadband 
once a motor is fired, i.e. 

(2-21) dg^^ = dCl-H) (degrees) 

Phase lags and delays in the system are treated as simple time delays. 
They include control electronics or computer lags, gyro lag, control motor 
thrust response and decay times and an incremental effect due to bending 
cycles. These time delays are applied to the control motor thrust tum-on and 
turn-off after sensing the crossing of a switching point defined by equations 
(2-20) and (2-2l). Bate terms have the most significant effect on the error 
signal dynamics. The rate gyros are assumed to have a time delay based on a 
second order response. 

(2-22) tRG = ^RG/'^fRG 

where, 4 rg is ihe gyro damping factor 

fpQ is the natural frequency (hertz) 

An effective time delay induced by large amplitude bending mode oscillations 
following ignition is, 

(2-23) = ~ 2f” (seconds) 

where, rj is an arbitrary phase angle at the switching point 
f^is the bending frequency (hertz) 


- 0O -n t)dt 


I i= 


= 0 


t < 


^1 P K;^edt ; |K;^e'iil<GLIMIT 

1“ *^LIMIT SIGN ) > j 
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Bending can increase or decrease the time delay. However, in the routine if 
the bending effect reduces the total system delay to a time advance, the time 
delay is set to zero. Total control motor thrust tum-on and turn-off delay 
time is, 


(2-24) 


^on = + "tHG + % + 

■toff = tG + tRG + + toffjjj 


Cseconds) 

(seconds) 


where. 


tQ is the guidance electronics delay (sec) 

tgj, is the control motor effective turn-on delay (sec) 

^offm control motor effective turn-off delay (sec) 


2.3 Yaw Axis Equations 


¥hen analyzing capture maneuvers in the yaw plane the equations used 
are basically the same as for pitch. However, horizontal winds from the side 
are not reduced by the sine of the flight path angle as in pitch. By entering 
the flight path angle as 90 degrees in input the correct aerodynamic moments 
will be generated for a sidewind. Equation (2-13) degenerates to. 


(2-25) = 57.3 V / Vv (degrees) 

where, a^, is angle of attack in yaw and is sidewind velocity. 

Plight path error becomes a flight path error measured in the yaw plane. 

Control motor designations become, 

(2-26 ) Ft right yaw motor 

Fr left yaw motor 

Sign convention trades the pitch-nose-up positive direction for a yaw 
nose-right positive direction. 


2.4 Roll Axis Equations 

Roll axis equations of motion can be simulated by appropriate choice of 
input variables. Sign convention is shown in Figure 3* Roll accelerations 
are, 

(2-27) ^ (deg/sec2) 

Roll control acceleration is, 

(2-28) = 57.3(Ft - Fb)Zc/ 12 (deg/sec2) 


- 7 - 



where, Fj ia roll right control motor thrust (lbs) 

Fg is roll left control motor thrust (lbs) 

Zq is moment arm about centerline (inche?) 

Ij is roll moment of inertia (slug-ft^) 

This ia similar to Equation (2-2) if the control motor cant angle (A^) is 
input as 90 degrees* 

Aerodynamic disturbances due to a nonlinear rolling moment coefficient 

iSf 

Cdeg/sec2) 


Cft 2/deg) 


(2-29) =57.3 CiSdrefq/12 


where. 


which is of the form of equations (2-3) and (2-5)* The following input data 
must correspond to pitch channel variables as follows, 


(2-31) 


Pitch 

(■CNaS)a=o 

Roll e 
/5CiS] 

\ S<^ 4>‘ 

*cg 

0 

^cp 

-dj,gf 

y 

0, 

•Vw 

0. 

V 

1. 


Booster induced roll torque accelerations are, 

(2-32) 57*3 I* / (deg/sec2) 

where L is the roll torque in (ft-lbs). This is similar to equation (2-4) 
with the following input changes. 

Pitch Roll 

T . 1 . 

fr L (ft-lbs) 

X 687.6 

0 . 

The changes in input data for the roll axis are defined in paragraph 
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2.5 Statistical Analysis 


Treatment of the random variables for input and output is discussed in 
the • following paragraphs. 


2.5.1 Input Handom Variables 


Random variables used in the analysis include mass properties, 
trajectory deviations, control system characteristics, thrust misalignment and 
wind. Most of these variables are treated as Gaussian (Normal. Distribution 
Function) and are entered as a mean and standard deviation. These include: 


moment of Inertia 
center of gravity station 
booster thrust 
dynamic pressure 

control system switching hysteresis 

control motor tum-on and turn-off delay times 

rate to displacement gain ratio 

control motor thrust level 

guidance electronics delay time 

rate gyro natural frequency and dampii^g 

initial rate 

wind speed 

flight path angle 

booster thrust misalignment 


Random samples of normally distributed parameters are computed as: 

(2-53) yi =• + RND a 

where (RND) is the Random Normal Deviate. 


Initial attitude error is computed based on the previous stage 
sensitivities to angle of attack and aerodynamic dissymmetries with respect to 
angle of attack. An independent contribution to initial attitude error is 
input as a Log-Normal Distribution of an oscillation amplitude. From Equation 

2-15, 

= ^amp ^in rj 

The amplitude of the sine wave is a positive number having a minimum value of 
zero. Several distribution functions can be used to describe the amplitude 
(i.e., Log-Normal, Rayleigh, Poisson). The Log-Normal was assumed since it is 
easily handled in the computer program by the Normal Distribution 
pseudo-random number generator. To accomplish this a new variable, *y'» is 
defined in input, 

y = ^amp 
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The mean and standard deviation of (y) the natural logarithm of amplitude is 
input * In the routine (y) is treated as Gaussian* The amplitude is then 
computed by 

(2-34) ^amp = + RND a ) 

where, fi is the mean of y 

a is the standard deviation of y 
END is the Random Normal Deviate 

The incremental attitude error at ignition includes a random phase angle (17) . 
Internally this phase is assumed to be a Uniform Distribution function ranging 
from - rr/2 to /2 radians. It is computed as 

( 2 - 35 ) 1 - ^ ( Q - 0.5) ; 0<Q<1 

where, Q, is a pseudo-random number ranging from zero to one. 


2.5.2 Output Statistical Analysis 


Individual capture transient information is stored in arrays. These 
include maximum overshoot attitude error, time to reach this maximum , and the 
total impulse expended by the control motors. These arrays are arranged in 
ascending order. Discrete probability levels are computed by Equation (2-36) 
and output with each sample. 


( 2 - 36 ) 

where y^^ 

% 

N 



is the value of the parameter 

is the location in the array of y 

is the total number of samples computed 


Histogram data are also output in tenas of the number of values falling 
in each cell specified in the input. Cell location algorithms used are; 

(2-37) HIST (ICELL) =* HIST(ICELL) + 1 

where, HIST (ICELL) contains the cumulative number of values in the 'i'th cell 


and, 

(2-38) ICELL = (yi - ymin)/Ay + 1 

y^ is the value of the parameter being tested 
ymin lowest range of the histogram specified in input 

(usually zero) 

Ay is the cell width computed from the input histogram range 
and number of cells. 

Capture parameters computed are always positive numbers. Therefore the 
mi nimuTn expected value (ymin) is zero. If a value greater than the maximum 
specified histogram is computed it is put into the last cell. 
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A divergent case will not te included in the discrete probability- 
tables. The number of such cases are printed out following the discrete 
probability tables. If there are any divergent cases the probability of 
successful capture can be computed externally. 

(2-39) P(SUCCESS) = l-P(FAILUEE) 

(2-40) P( FAILURE) = Np/N^QT 

where Np is the number of failures 

NipQiji is the total number of samples 


If there are capture manuevers which are considered failures due to 
high overshoot angles the probability of failure is calculated from Equation 
(2-40) where 


(2-41) Np = Np/c Nd 

(2-42) %0T = N + Njj 


where Np^g 


Nd 

N 


are the number of cases which captured but are considered 
failures 

is the nimber of divergent cases calculated 

is the sample size requested (note this does not include 

divergent cases or other error modes encountered) 


Generally a finite sample size computed by Monte Carlo techniques makes 
application of Equation (2-40) somewhat inaccurate. Therefore it is 
desireable to attach a confidence level to the probability. This can be done 
by using the Binomial Failure Rate Tables. These are based on the solution of 
the following equation for confidence coefficient: 


(2-43) 


Nf 

/3=i-y 


i=o 


Ntot! 

j • ( Ntot” J ^ • 



A |*T0T" j 
'■) 


where )8 is the confidence coefficient 

Np is the total number of failures 

Nijiqj is the total number of trials 

F is the upper confidence limit on the true probability of 
failure 


A table of the upper confidence limit of the true probabilty of failure is 
presented in Figure 5 for sample sizes of 200, 500 and 1000. 
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3.0 PROGRAM DESCRIPTION 


3.1 General 


This computer program is coded in FORTRAN IV for a CDC CYBER 175 
system. Code is compatible with ANSI standards. It ia arranged to operate 
with standard card input and line printer output. 

A main routine and four subroutines require approximately 8K words of 
computer memoiy. Program flow and user instructions are presented in the 
following paragraphs. Input and output of sample problems are illustrated 
along with the detailed descriptions. These can be used to checkout the 
program on another machine. 


3.2 Progreim Plow 


Program flow is straight forward in fiye basic parts, 

. input data 
. initilize constants 
. compute Monte Carlo capture samples 
. compute statistical distribution 
. output Monte Carlo results 

All four subroutines are single level. Blank common is used to pass 
information between MAIN and subroutines RNDX, and PAGEHD. 

A flow chart showing the general flow of the MAIN Program is presented 
in Figure 6. More detailed charts of specific branches of the routine are 
presented in Figure 7. A complete list of the FORTRAN source program is 
presented in the Appendix. 

A description of the subroutines is presented in the following 
paragraphs. 


3.3 Subroutine Description 

Four subroutines are used to support the MAIN program; ASCEND, PAGEHD, 
RNDX, and TBLU. A brief description of each is presented below. 

ASCEND 


This subroutine rearranges an array of numbers in monotonically 
increasing or decreasing order. The call to this subroutine is: 
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CALL ASCEND (L, VAL, M) 


L - is the number of values in array VAL (input) 

VAL - is the array to be rearranged (input and output) 

M - is input option 

M = 0 for ascending order 
M = 1 for descending order 

PAGEHD 

This subroutine ejects a page, prints run number page number and title 
infonnation at the top of the page. 

The call statement is 


CALL PAGEHD 

Page number, line number, run number and title information are 
transferred by blank common from the MAIN routine. 

RNDX 

This is a pseudo random number generator for a Uniform or a Normal 
distribution function. It uses a seed integer 'K' to compute a uniformly 
distributed number 'Q' which ranges from zero to one. If the Normal 
Distribution is required the Uniform Distribution U(0,l) is transformed to the 
Normal Distribution N(0, l). A random normal deviate, 'RXD', is obtained from 
the value of 'Q'. The call to this subroutine is, 

CALL RNDX(l) 

L = 0 gives random normal deviate N(0, l) 

L = 1 gives uniform U(0, l) 

The values of 'Q', 'RXD' and the changed random number seed value 'K' 
are located in blank common. 


TBLU 

This is a single table lookup subroutine. It is based on linear 

interpolation between points from a single array having alternating values of 
abscissas and ordinates. Abscissas must be in ascending order. If the 
abscissa given is out of the range, either the first or the last table 
ordinate is returned depending on the direction of the abscissa. 

The call statement is 

-CALL TBLU (NT, Y, X, T, M) 

NT - number of values in table 'T' including abscissas and ordinates 
(input) 

Y - is the ordinate found (output) 

X - is the given abscissa (input) 

T - is the table of alternating abscissas and ordinates (input) 

M - is the table locator index. 'M' must be greater than zero and 

less than 'NT' . 'M' returns the index for the current location 

found for the abscissa and should be used for the next table 
lookup of the same table to minimize search time. 
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3.4 Error Modes 


There are several error Modes which can be encountered during the 
numerical integration of a capture msinuever. These result in cases which 
cannot be completed within the input limitations or the computational 
procedure. A description, of each of these is given below along with the test 
equations. Printout of some key random variables is made when an error mode 
is encountered. 


3.4.1 Deadband Overahoot (Error Mode l) 


This error occurs when the error signal over a single integration step 
crosses a tum-on and turn-off point without turning off a control motor (see 
Figure 8) . This is a result of either a poorly designed control system or too 
coarse of an integration step size. If this is encountered it is suggested 
that integration step size be reduced significantly. If this does not affect 
the results, recheck all input data. Improper input data is usually the 
cause. As an example if the bending frequency is arbitrarily set to zero or 
some other low value the effective time delay will become very large. 

Incorrect units on input data is another common error. 

If no errors cein be found, the control system, disturbing 
accelerations, control gains, and deadbands should be reviewed along with the 
specified system time delays. Reduction of control force, widening deadbands 
or reducing gains should reduce the number of deadband overshoot cases. 


3.4.2 Divergence 


Divergence is the loss of attitude control due to excessive disturbance 
(see Figure l) . It occurs when the attitude error reaches the limiting value 
determined by input data. These cases are not included in the end statistical 
results and must be treated as described in paragraph 2.3*2. 


3.4.3 Time Limit 


In order to restrict computer run time the user must input a time limit 
for the capture maneuver. If this were not done certain conditions could lead 
to an endless or excessively long capture maneuver which in most cases would 
be trivial. As an example if there were very low disturbances and rate it 
would take a relatively long time to reach a deadband to fire a control 
motor. By inputing a realistic time limit such as 5 or 10 seconds, these 
trivial cases would be terminated and an Error mode would be encountered. 
Review of output data generally would reveal the seriousness of the error. If 
an excessive number of these errors are encountered either the input data is 
in error or the time limit specified is too short for the control system. 


- 14 - 



3*5 Input Data Description 


Input data descriptions for pitch, yaw, or roll analyses are described 
in the following paragraphs. Pitch, yaw and roll options use the same 
routine. Only the input data is redefined. Sample problem input listings are 
presented in Figures 9 through 11 to aid in the descriptions. Input data 
occurs in eight (S) basic groups: 

(1) card of integers which includes option designations, sample size 
and run number, 

(2) two cards of arbitrary hollerith data which is printed out for 
identification, 

(3) a card containing an integer initial random number generator seed, 

(4) nineteen (l9) cards containing the floating point values of the 
mean and standard deviation of random constants, 

(5) nineteen (l9) cards containing constants which include limits, 
integration data and non-random constants, 

(6) three cards for integral gain option if used, 

(7) four (4) tables of time histories, 

(8) optional histogram output information. 

Groups (1), (2), (3), (6), and (8) are the same for pitch, yaw, or roll 
axes. Groups (4), (5)» and (7) are generally the same for a pitch or yaw axis 
analysis with minor differences in definition. Roll axis analysis results in 
many definition changes which are described in paragraph 3»5.9» 


3.5.1 First Card (Group l) 


The first input card contains eight fixed point numbers entered in 
fields of five. These must be entered right justified. Definitions of each 
of these constants is given below. 

Card Column FORTRAN Variable Description 

1-5 NRUN An arbitrary run number 

10 HPRINT Output print option 


0 print out summary data 

1 also print values at each 
overshoot on each case 

2 Same as (l) with addition of 
thrust on and off switchiiag 
points 

3 print out data at each 
integration step of each case 

11 - 15 NTRUN Number of random samples to be 

computed by the Monte Carlo 
technique. (Current dimensions 
limit this to 1000 samples) 
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FORTRAN Variable 


Deacription 


INC 


NSTAGE 


N 


NHIST 


NKI 


Incremental number of samples 
for computing the running mean 
and standard deviation Of the 
capture results 

Desired control system state at 
start of the capture 

2 Control system activation 
occurs simultaneously with 
ignition or separation 
event. Control motors are 
off initially even if error 
is outside of deadband 

3 Control system activated 
prior to ignition. If 
initial conditions are 
outside of deadband control 
motor thrust is applied 
immediately 

Number of overshoots desired 
before the capture maneuver is 
assumed to be completed. If an 
integral gain is used (NKI=l) , 
overshoots are counted only 
after the integral gain is 
switched into the loop 

Option to obtain histogram data 
for capture parameters 

0 No histogram data, capture 
parameters are printed out in 
unarranged order of 
computation 

1 Histogram data is output, 
(Group 8 input is required) 

Option for integral of attitude 
error gain function in control 
logic 

0 No integral gain term 

1 Integral gain used (Group 6 
data must be input) 



3. 5* 2 Arbitrary Identification (Group 2) 


The second and third cards contain user defined arbitraiy labeling 
information. The alphanumeric information is contained on the first 
seventy-eight (78) columns of each card. It is output at the top of the page 
in two lines the same as input. Internally the information is stored in array 
NTITLE which contains 26 words each containing 6 characters. 


3 . 5.3 Initial Random Number Seed (Group 3) 


The fourth card contains the initial value of the random number seed 
integer (K). This should be limited to a seven (?) digit integer located in 
columns 1 through 10 right justified. 


3 . 5*4 Random Constants (Group 4 ) 


Cards (5) through (23) contain the mean (XM) and standard deviations 
(XSD) of nineteen random constants. These are read in two fields of fifteen 


(FORMAT 2E15.8). These 

are: 



Card No. 

FORTRAN 

Variable 

Symbol 

Units 

Description 

5 

XM(1), XSD(l) 

ly 

slug-ft^ 

Moment of inertia, mean 
and standard deviation 

6 

XM(2), XSD(2) 

^cg 

inches 

Center of mass body 
station location, mean 
and standard deviation 

7 

3ffl(3), XSD(3) 

^/"^nom 


ratio of booster motor 
thrust to nominal value, 
mean and standard 
deviation 

8 

XM(4), XSD(4) 

9./ ^nom 

— 

ratio of dynamic pressure 
to nominal value, mean 
and standard deviation 

9 

XM(5), XSD(5) 

d 

degrees 

deadband halfwidth, mean 
and standard deviation 

10 

XM(6), XSD(6) 

H 

— 

switching hysteresis 
(l"<ioff/'ion) 

11 

XM(7), XSD(7) 

^onm 

seconds 

control motor turn-on 
delay time 

12 

XM(8), XSD(8) 

^offn, 

seconds 

control motor turn-off 
delay time 



Card No. 

FORTRAN 

Variable 

Symbol 

Units 

Description 

15 

XM(9), XSD(9) 

Kr/Kd 

seconds 

rate to displacement gain 
ratio 

14 

XM(lO), XSD(IO) 

Fij 

lbs 

upper control motor 
thrust (right control 
motor for yaw) 

15 

XM(ll), XSD(ll) 

Fb 

lbs 

lower control motor 
thrust (left control 
motor for yaw) 

16 

XM(12), XSD(12) 


seconds 

control electronics time 
delay 

17 

XM(13), XSD(13) 

^RG 

— 

rate gyro damping factor 

18 

XM(14), XSD(U) 

^‘rg 

hertz 

rate gyro natural 
frequency 

19 

XM(15), XSD(15) 


degree 

natural log of attitude 
error oscillation 
amplitude prior to 
ignition 

20 

XM(16), XSD(16) 

^0 

deg/sec 

attitude rate initial 
condition 

21 

XM(17), XSD(17) 

Vw 

ft/sec 

headwind velocity at 
ignition (right sidewind 
for yaw) from wind 
profile at ignition 
altitude 

22 

XM(18), XSD(18) 

Vw 

ft/sec 

headwind velocity at 
second time point, card 
35 in Group 5 (right 
sidewind for yaw) from 
wind profile at 
appropriate altitude 

23 

XM(19), XSD(19) 

^0 

degrees 

flight path angle at 
ignition measured from 
local horizontal (yaw 
path angle usually 90 
degrees mean) 
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3«5«5 Constants (Group 5) 


This group of data includes non-random constants for the integration, 
control system characteristics, vehicle and trajectory characteristics. The 
first two cards in this grouj) contain aerodynamic coefficients, two to a card 
with Format (2E15.8). These are followed by seventeen (l7) cards each having 
a single value read in with Format (E15*8). They are: 


FORTRAN 


Card No. 

Variable 

Symbol 

Units 

Description 

24 

CNAS 

(CNoS)a=0 

ft^/deg 

aerodynamic normal force 
coefficient slope at a=0; 


A 3 

a3 

ft^/deg5 

cubic term in aerodynamic 
normal force coefficient 

25 

XCPO 

^cpo 

inches 

body station of 
aerodynamic center at a=0 


DCPDA 

5X(jp/Sa 

in/deg 

rate of change in 
aerodynamic center with 
angle of attack 

26 

XF 

Xc 

inches 

body station of control 
force 

27 

ZC 

Zc 

inches 

control motor location 
off centerline 

28 

XL AM 

^c 

degrees 

forward cant angle of 
control motor thrusts 

29 

XN 

Xr 

inches 

body station of booster 
nozzle throat 

30 

WB 

fb 

hertz 

bending frequency (use 
very large number if 
rigid body, e.g., lOOOOO) 

31 

TI 

to 

seconds 

start time (usually zero) 

32 

OMAX 


seconds 

present cutoff timer for 
capture integration. If 
this time is reached an 
error mode is indicated 
and data for this case is 
printed out 

33 

DELTO 

— 

seconds 

coarse integration step 
size 

34 

DELTMN 

— 

seconds 

minimum integration step 
size used when isolating 


switching points 
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FORTRAN 


Card No 

. Variable 

Symbol 

Units 

Description 

35 

AT 


seconds 

second table time used 
for velocity and wind 
velocity input 

36 

VO 

Vo 

ft/sec 

vehicle velocity at 
ignition 

37 

VI 

Vl 

ft/sec 

vehicle velocity at 
second time point (card 
35) 

38 

GDOT 

n 

deg/sec 

pitch torqueing rate or 
flight path angle rate 
(usually zero for yaw) 

39 

ASEN 

Ko 


sensitivity of attitude 
error to angle of attack 
for previous stage trim 
conditions 

40 

THCMO 

^e bias 

degrees 

attitude error bias of 
previous stage at zero 
angle of attack 

41 

TLIMIT 

% limit 

degrees 

maximum allowable capture 
overshoot angle. If this 
limit is reached the case 
is terminated and a 
divergence error is 
indicated 

42 

ADES 

®nom 

degrees 

predicted angle of attack 
on an undisturbed 
trajectory 


3.5.6 

Integral 

Gain Values 

(Group 6) 

Three (3) cards of constants are needed when 
control logic is used. These are input only if 
or greater. One value per card is input with a 

an integral of attitude error 
NKI on the first card is one 
FORMAT (E15.8) as follows: 

Card No. 

FORTRAN 

Variable 

Symbol 

Units 

Description 

43 

GINT 

Kl 

l/sec 

integral of attitude 
error gain 

44 

XLIM 

® limit 

degrees 

integral gain term 
limiting 

45 

TKI 

*KI 

seconds 

time after ignition at 
which integral gain term 
enters control logic 
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3.5*7 Time History Tables (Group 7) 


Four (4) tables of parameter time histories are input. These include 
dynamic pressure, booster thrust, the mean and standard deviation of booster 
thrust misalignment angle. Each of these tables is input with FORMAT 
(I5/(6E10.3) ) • The first card of each table contains the integer number of 
values in the array; it must be located in columns one (l) through five (5) 
and be right justified. Succeeding cards contain alternating values of time 
and ordinate with six (6) numbers per card in fields of ten (lO). Each of 
these tables are dimensioned to carry up to fifty values or twenty-five (25) 
time points each. 


FORTRAN 
Table Name 

Ordinate 

Symbol 

Units 

Absciaas 

Symbol 

Units 

Description 

TABQ 

'Inom 

lbs/ft2 

t 

sec 

nominal dynamic 
pressure 

TABT 

T 

lbs 

t 

sec 

booster thrust 

ETMEAN 

(r 

degrees 

t 

sec 

mean value of 
booster thrust 
misalignment 

ETSIG 


degrees 

t 

sec 

standard deviation 
of booster thrust 
misalignment 


3.5*8 Histogram Information (Group 8) 


This last group of input is used only when NHIST is one (l) or 
greater on the first card. This results in histogram cell counts for the 
overshoot angle, time to overshoot, and impulse expended output parameters. 
Input includes two cards of data. The first card contains the integer number 
of cells for the histogram (NCEL). This is entered in the first five (5) 
columns and is right justified. A maximum of 150 cells is included in the 
dimension statement. 

The second and last card contains the minimum and maximum values 
which specify the ranges of each histogram. These are entered with format 
(6E10.3) and are: 
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Card FORTRAN 


Column. 


Variable 

Units 


Description 

1-10 


CMNH 

degrees 


attitude error overshoot 
angle lower range (usually 
zero) 

11 - 20 


CMAXK 

degrees 


attitude error overshoot 
upper end of range 

21 - 30 


THINK 

seconds 


capture time lower end of 
range (usually zero) 

I-* 

1 

O 


TMAXH 

seconds 


capture time upper end of 
range 

41 - 50 


PMINH 

lb-sec 


capture impulse lower end of 
range (usually zero) 

51 - 60 


PMAXH 

lb-sec 


capture impulse upper end of 
remge 

In output if any data ia 
included in the first and last 

out of the range 
cells. 

specified they will be 



3 . 5.9 Roll 

Axis Input Data Differences 

When using the routine for roll azes analyses the definition of some 
terms is different. A sample problem input is presented in Figure 11. The 
changes from a pitch axis run are reflected in groups 4, 5» and 7. Since wind 
and angle of attack are not included in the roll azis equations certain 
parameters associated with these must be entered as a specified constant. 

These are denoted by asterisks in the description column in Figure 11. Group 
4 and 5 data are the same as presented in paragraphs 3*5.4 and 3*5.5 with the 
following exceptions. Asterisks indicate specified numbers must be entered. 

Card No. 


FORTRAN 

Variable 

Symbol 

Units 

Description or 
Specified Numbers 

6 


XM(2), XSD(2) 

*** 

•*« 

0. ^ 0 . 

7 


XM(3), XSD(3) 

*«* 

««« 

1., 0. 

14 


XM(IO), XSD(lO) 


lbs 

roll right control motor 
thrust. If a couple 
this must reflect both 
motors 

15 


XM(ll), XSD(ll) 

Fb 

lbs 

roll left control motor 
thrust 

21 


XH(17), XSD(17) 

»** 

*«» 

0., 0. 

22 


XM(18>, XSD(18) 

*«* 

*** 

0. , 0. 

23 


XM(19), XSD(19) 

«*« 

»»» 

0., 0. 
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Card No. 

FORTRAN 

Variable 

Symbol 

Units 

Description or 
Specified Numbers 

24 

CNAS 


ft 2/deg 

aerod 3 niamic rolling 
moment coefficient slope 
at roll error of zero 


A3 

^3 

ft2/deg^ 

cubic term in rolling 
moment coefficient 

25 

XCPO, 

~^ref 

inches 

negative of the 
aerodynamic reference 
length 

26 

XF 

*** 

«*« 

0. 

27 

ZC 

Zc 

inches 

roll control motor 
location from vehicle 
centerline 

28 

XL AM 

*«* 

*** 

90. 

29 

XN 

«*• 

’WWW 

687.6 

30 

WB 

fb 

hertz 

torsional mode frequency 

35 

AT 

*** 

««* 

1. 

36 

VO 

*** 

«*« 

1. 

37 

VI 

*** 

«»« 

1. 

38 

GDOT 

a 

deg/sec 

commanded roll attitude 
rate 

39 

ASEN 

«** 

•«« 

0. 

40 

THCHO 

0e bias 

deg 

initial roll attitude 
error bias 

42 

ALES 

««« 

««« 

0. 


Group 7 tables are redefined from those presented in paragraph 3*5»7* 
The dynamic pressure is the same. The inputs are, 


FORTRAN 

Ordinate 


Absciaas 



Table Name 

Symbol 

Units 

Symbol 

Units 

Description 

TABQ 

9nom 

lbs/ft2 

t 

sec 

nominal dynamic 
pressure 

TABT 

*«* 

1. 

t 

sec 

unity multiplier 
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FORTRAN 
Table Name 

Ordinate 

Symbol 

Units 

Absciaas 

Symbol 

Units 

Description 

EMEAN 

L 

ft-lbs 

t 

sec 

mean value of 
booster induced 
roll torque 

ETSIG 

L 

ft-lbs 

t 

sec 

standard deviation 
of booster roll 
torque 



3.6 Output 

Data Description 



All output for this program is printed on output tape unit six (6). 

The amount of output is controlled by options selected in the first card of 
input (paragraph Output is made in groups, i. e., 

1) running mean and standard deviation and individual sample data 

2) unarranged capture overshoot, time and impulse for all cases 
(NHIST=0 only) 

3) histogram cell data for capture overshoot, time and impulse 
(NHIST =1) 

4) discrete probability versus overshoot, time and impulse arranged in 
ascending order 

5) summary of divergent cases, abnormal terminations and final pseudo 
random number seed 

Detailed descriptions of the output with sample problems are presented 
in the following paragraphs. 


3.6.1 Running Mean, Standard Deviation and Sample Data (Group l) 


The amount of data output in this first section depends upon the 
selected options. The amount of information output is determined by the 
parameter NPRINT on the first card of input described in paragraph 3* 5.1. 

NPRINT = 0 

When NPRINT = 0 the running mean and standard deviation is printed 
after the specified increment (INC) or random captures. The sample output for 
this option is shown in Figure 12 for a 199 sample run with INC = 50. The 
following nomenclature applies. 
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ITER NO 

is the current sample at which all previous normally 
computed cases are included in the mean (MEAN) and 
standard deviation (SlfflA) calculations. 

THETA 

is the absolute value of the maximum attitude error 
overshoot in degrees. 

TIME 

is the time to reach the maximum overshoot angle 
(seconds) , 

IMPULSE 

is the control motor total impulse expended during 
capture in pound-seconds 

Also interspersed within these output lines will be four lines 
indicating abnormal terminations as described in section 3*4. The sample 
problem presented in Figure 12 had one divergent case between sample 150 and 
199. Output for the abnormal case is always as shown in Figure 12 with the 
following difinitions; all values occurred at time of termination. 

TIME (SEC) 

is the time from stage ignition at which the error 
occurred 

RATE (DEG/SEC) 

is the attitude angular rate 

ATT ERROR (DEG) 

is the attitude error 

Q (PSF) 

is the dynamic pressure 

ALPHA (DEGREE) 

is the angle of attack or roll attitude error for a 
roll capture 

THRUST MIS (DEG) 

is the booster thrust misalignment in degrees or the 
roll torque in ft-lbs 

F TOP (LB) 

is the upper pitch, right yaw, or roll right control 
motor thrust 

F LOVIER (LB) 

is the lower pitch, left yaw, or roll left control 
motor thrust 

LAST OVERSHOOT 

la the maximum overshoot angle (zero rate crossing) 
prior to the termination 


Above information is output for any printout option and is interspersed 
with other outputs in Group 1. 


NPRINT = 1 

This selection of output includes all overshoots which occurred during 
the capture maneuver. A sample output is shown in Figure 13. Three 
overshoots were selected for this run (N parameter on first card of input). 
Definitions of terms include, 

THETA is attitude error in degrees 

THETADOT is attitude rate in degrees per second 
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TIME 


is time from stage ignition in seconds 

IMPUI5E is control motor impulse expended in pound-seconds 

NPRINT =• 2 

With this option the Group 1 output includes all items discussed above 
plus points in each capture at control deadband crossing and thrust-on and 
thrust-off points. A sample output is shown in Figure 14. 

NPRINT = 3 

This option increases the amount of output for each sample to include 
all integration steps. A sample output is presented in Figure 15. It is 
helpful for individual capture maneuvers. As shown on the second and third 
page of Figure 15 one sample ended in divergence. 


3.6.2 Unarranged Parameters (Group 2 - NHIST = O) 


In ceirtain circumstances, the order in which attitude error overshoot, 
capture time, and impulse are desired as well as the relationship of these 
three parameters. This can be obtained by setting NHIST = 0 on the first card 
of input. It prevents the histogram output from being generated and output 
also. A sample output is shown in Figure 16. 


3.6.3 Histogram' Cell Data (Group 3 - NHIST ■» l) 


Capture overshoot angle, time, Eind impulse are counted and located in 
histogram cells if NHIST = 1 on the first card of input. The number of cells 
and ranges for the histogram are user selected in the Group (8) input 
described in paragraph 3*5.8. A sample output is presented in Figure 17. The 
following is output for each cell. 


CAPTURE ANGLE the attitude error overshoot angle range of 

the cell in degrees 


NO the number of cases which had values in this 

cell range 

TIME the time to maximum overshoot cell range in 

seconds 


IMPULSE 


the control motor impulse expended in 
pound-seconds 
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5.6.4 Discrete Probability (Group 4) 


All options include this output. It includes the maximum overshoot 
angle, capture time, and impulse expended versus probability level as defined 
by Equation 3-1, 

(3-1) P = %/ (NTRUN + 1) 

where, P is the probability of not exceeding the given value 

is the sample number after being arranged in ascending 
order for the particular variable 

HTRUN is the number of normally completed cases, specified on the 
first card 

A sample problem output is presented in Figure 18. It should be noted 
that the probability of successful capture must include these results plus the 
abnoimally tenninated cases as discussed in paragraph 2.5.2. 


3.6.5 Summary of Abnormal Cases (Group 5) 


This page of output occurs with all options. A sample is presented in 
Figure 19. It includes the total number of divergent maneuvers, as well as 
those which encountered a "deadband overshoot" or a time limit as discussed in 
Section 3.4. The final integer value of the pseudo-random number generator 
seed is also output on this page. This number can be used as an input to a 
successive run to increase the sample size based on a pseudo-random number set 
which is "uncorrelated" from the previous run. 
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Figure 1 

Hiase Plane Trajectories of Capture Haneuver 
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Figure 2 

Control System Block Diagram 



q , a , y 















Figtore 3 
Sign Convention 


PITCH AXIS 



ROLL AXIS 
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Figure 4 

Atypical Aerodynamic Coefficients 

SCOUT Second Stage~Mach 4 
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Figure 5 

Tables of Binomial Upper Confidence Limits 
Probability of Failure - 200 Samples 


Confidence 


Coefficient .500 

,800 

,900 

.950 

.990 


.995 

No . of 

fai lures 
0 

0.003460 

0,008015 

0.011447 

0,014867 

0.022763 

0 

,026144 

1 

0.008377 

0,014897 

0,019309 

0,023498 

0,032728 

0 

,036559 

2 

0,013348 

0,021274 

0.026391 

0,031143 

0,041362 

0 

,045534 

3 

0*018330 

0,027402 

0.033098 

0.030310 

0*049348 

0 

,053803 

4 

0.023315 

0.033377 

0,039570 

0.045180 

0.056929 

0 

,061631 

5 

0,028303 

0,039243 

0,045879 

0,051843 

0.064224 

0 

.069146 

6 

0,033292 

0,045028 

0,052064 

0,058350 

0.071304 

0 

,076425 

7 

0,038282 

0,050749 

0,058153 

0,064733 

0,078212 

0 

,083516 

8 

0,043272 

0.056417 

0,064162 

0.071014 

0,084981 

0 

.090455 

9 

0.048263 

0,062041 

0,070104 

0,077211 

0,091632 

0 

,097264 

10 

0.053254 

0,067620 

0,075990 

0,083335 

0,098182 

0 

.103960 

11 

0,058245 

0.073102 

0.001826 

0.089396 

0.104640 

0 

.110560 

12 

0.063235 

0,070707 

0.087618 

0,095401 

0.111030 

0 

.117080 

13 

0.068227 

0,084207 

0,093371 

0,101360 

0.117340 

0 

.123520 

14 

0,073218 

0.009683 

0,099089 

0,107270 

0,123600 

0 

.129890 

15 

0,078209 

0.095139 

0*104780 

0,113140 

0.129800 

0 

.136200 

16 

0,083200 

0,100580 

0,110430 

0,118970 

0.135940 

0 

,142460 

17 

0,088192 

0,105990 

0,116060 

0.124770 

0.142040 

0 

.148660 

18 

0,093183 

0,111400 

0,121670 

0,130540 

0,148090 

0 

,154810 

19. 

0,098174 

0,116780 

0,127250 

0,136280 

0.154110 

0 

,160920 

20 

0.103170 

0,122160 

0,132820 

0,141990 

0.160080 

0 

,166980 

21 

0,108160 

0,127520 

0,138360 

0,147680 

0,166020 

0 

.173010 

22 

0,113150- 

0.132870 

0,143880 

0,153340 

0.171930 

0 

,179000 

23 

0,118140 

0.138200 

0,149390 

0.158980 

0.177800 

0 

.184950 

24 

0.123130 

0,143530 

0,154880 

0.164590 

0,183640 

0 

.190880 

25 

0,128120 

0,148840 

0.160350 

0.170190 

0.189460 

0 

,196770 

26 

0,133110 

0,154150 

0.165810 

0.175770 

0.195250 

0 

,202630 

27 

0,138110 

0,159440 

0.171250 

0.181330 

0,201010 

0 

,208460 

28 

0,143100 

0,164730 

0. 176680 

0*186870 

0.206750 

0 

,214270 

29 

0,148090 

0.170010 

0,182100 

0,192400 

0,212470 

0 

.220050 

30 

0,153080 

0,175280 

0,187500 

0,197910 

0.218160 

0 

,225800 

31 

0,158070 

0,180540 

0,192900 

0,203410 

0,223830 

0 

.231530 

32 

0,163060 

0,185790 

0,198280 

0,208890 

0,229480 

0 

.237240 

33 

0,168060 

0,191040 

0,203650 

0.214350 

0,235110 

0 

,242930 

34 

0,173050 

0. 196280 

0,209010 

0.219800 

0,240730 

0 

.248590 

35 

0,178040 

0,201510 

0,214350 

0.225240 

0,246320 

0 

,254240 

36 

0,183030 

0,206740 

0,219690 

0.230670 

0.251890 

0 

.259860 

37 

0,188020 

0,211960 

0.225020 

0,236080 

0.257450 

0 

,265470 

38 

0,193010 

0,217170 

0,230340 

0,241480 

0.262990 

0 

,271060 

39 

0,198010 

0,222380 

0,235650 

0,246870 

0.268520 

0 

.276620 

40 

0*203000 

0,227580 

0.240950 

0.252250 

0.274020 

0 

,282180 

41 

0*207990 

0,232780 

0,246250 

0,257620 

0.279520 

0 

.287710 

42 

0.212980 

0,237970 

0,251530 

0,262970 

0,285000 

0 

.293230 

43 

0,217970 

0,243150 

0.256810 

0,268320 

0,290460 

0 

.298730 

44 

0.222960 

0,240330 

0.262070 

0.273650 

0.295910 

0 

,304210 

45 

0,227960 

0,253510 

0,267330 

0,278980 

0.301340 

0 

,309680 

46 

0. 232950 

0*258680 

0,272590 

0,284300 

0,306760 

0 

.315140 

47 

0*237940 

0,263840 

0.277830 

0.289600 

0.312170 

0 

.320580 

48 

0,242930 

0,269000 

0,283070 

0,294900 

0,317570 

0 

,326000 

49 

0,247920 

0,274160 

0.288300 

0,300190 

0.322950 

0 

.331420 

50 

0,252910 

0,279310 

0,293530 

0,305470 

0,328320 

0 

,336810 




Figure 5 (cont'd) 

Tables of Binomial Upper Confidence Limits 
Probability of Failure - 500 Samples 


Confidence 
Coefficient *500 

4800 

.900 

.950 

.990 


.995 

No. of 

failures 

0 

0*001385- 

0*003214 

0*004595 

0*005974 

0*009168 

0 

,010541 

1 

0*003354 

0*005977 

0*007757 

0.009452 

0*013202 

0 

,014765 

2 

0*005344 

0*008539 

0*010609 

0*012538 

0*016705 

0 

,018413 

3 

0.007339 

0*011002 

0*013312 

0*015434 

0.019949 

0 

.021781 

4 

0.009336 

0.013405 

0*015923 

0,018213 

0*023033 

0 

*024973 

5 

0*011333 

0.015766 

0*018470 

0*020910 

0*026005 

0 

,028042 

6 

O.OI 333 O 

0*018095 

0*020969 

0*023547 

0,028892 

0 

,031018 

7 

0*015328 

0*020399 

0*023430 

0.026135 

0,031713 

0 

,033921 

6 

0.017326 

0*022683 

0*025860 

0*028684 

0.034479 

0 

,036765 

9 

0*019325 

0*024949 

0*028264 

0,031200 

0.037200 

0 

.039558 

10 

0.021323 

0,027202 

0,030647 

0.033688 

0*039882 

0 

,042309 

11 

0.023321 

0*029442 

0*033010 

0.036153 

0.042530 

0 

,045023 

12 

0,025320 

6*031670 

0*035357 

0.038596 

0,045149 

0 

,047705 

13 

0*027318 

0*033890 

0*037690 

0,041020 

0,047741 

0 

*050358 

lA 

0*029317 

0*036100 

0*040009 

0,043427 

0,050310 

0 

*052984 

15 

0*031315 

0.038303 

0*042316 

0*045820 

0.052858 

0 

,055588 

16 

0*033314 

0*040499 

0,044613 

0*048198 

0*065387 

0 

,058170 

17 

0*035312 

6*042688 

0.046899 

0.050564 

0.057898 

0 

,060733 

18 

0*037311 

0,044871 

0*049177 

0,052918 

0*060392 

0 

,063278 

19 

0*039309 

0*047048 

0.051446 

0,055262 

0*062872 

0 

.065806 

20 

0*041308 

0*049221 

0.053707 

0.057596 

0*065338 

0 

.068319 

21 

0*043306 

0,051389 

0.055962 

0.059920 

0*067791 

0 

,070818 

22 

0*045305 

0,053552 

0.058209 

0*062236 

0*070232 

0 

.073303 

23 

0.047303 

0.055711 

0.060450 

0.064543 

0*072661 

0 

.075776 

24 

0.049302 

0*057866 

0.062685 

0.066843 

0.075080 

0 

.078237 

25 

0*051301 

0,060018 

0*064915 

0*069136 

0*077488 

0 

.080687 

26 

0*053299 

0*062166 

0.067139 

0*071422 

0*079887 

0 

,083126 

27 

0*055298 

0*064310 

0*069358 

0*073702 

0*062277 

0 

.085556 

28 

0*057297 

0*066452 

0.071572 

0*075975 

0*084658 

0 

,087975 

29 

0*059295 

0,068590 

0*073782 

0*078243 

0,087032 

0 

,090386 

30 

0*061294 

0*070726 

0.075988 

0.080505 

0*089397 

0 

,092789 

3-1 

0*063292 

0*072859 

0.078189 

0.082762 

0*091755 

0 

.095183 

32 

0*065291 

0*074989 

0.080386 

0.085013 

0*094106 

0 

,097569 

33 

0*067290 

0,077117 

0*082580 

0,087260 

0*096450 

0 

,099948 

34 

0*069288 

0*079242 

0*084770 

0,089502 

0*098788 

0 

.102320 

35 

0.071287 

0*081365 

0*086956 

0.091740 

0*101120 

0 

,104680 

36 

0,073286 

0*083486 

0,089139 

0*093974 

0.103440 

0 

.107040 

37 

0*075284 

0,085605 

0*091319 

0*096203 

0,105760 

0 

,109390 

38 

0.077283 

0,087722 

0.093496 

0.098428 

0,108080 

0 

.111740 

39 

0*079281 

0*089837 

0*095670 

0*100650 

0*110390 

0 

*114080 

40 

0*081280 

0*091950 

0*097841 

0,102870 

0*112690 

0 

*116410 

41 

0.083279 

0*094061 

0*100010 

0*105080 

0.114990 

0 

,118740 

42 

0*085277 

0*096170 

0*102170 

0*107290 

0*117280 

0 

,121060 

43 

0,087276 

0*098278 

0* 104340 

0*109500 

0*119570 

0 

,123380 

44 

0*089275 

0*100380 

0.106500 

0.111700 

0.121850 

0 

.125690 

45 

0,091273 

0* 102490 

0*108660 

0*113900 

0.124130 

0 

,128000 

46 

0,093272 

0,104590 

0,110810 

0,116100 

0*126410 

0 

,130300 

47 

0*095271 

0,106690 

0.112960 

0*110300 

0.128680 

0 

,132600 

48 

0*097269 

0,108790 

0,115110 

0*120490 

0*130950 

0 

.134890 

49 

0*099268 

0,110890 

0*117260 

0*122600 

0*133210 

0 

,137180 

50 

0*101270 

0*112990 

0.119410 

0.124070 

0*135470 

0 

,139470 

51 

0*103270 

0*115080 

0,121550 

0.127050 

0.137720 

0 

*141750 

52 

0.105260 

0,117180 

0,123700 

0.129230 

0,139980 

0 

,144030 

53 

0.107260 

0,119270 

0,125840 

0*131410 

0*142220 

0 

.146300 





Figure 5 (concluded) 

Tables of Binomial Upper Confidence Limits 
Probability of Failure - 1000 Samples 


Confidence 
Coefficient «500 

#800 

.900 


,950 


.990 

,995 

No. of 

failures 

0 

0 

000693 

0#001608 

0,002300 

0, 

002991 

0 

004595 

0,005284 

1 

0 

001678 

0.002991 

0.003884 

0. 

004735 

0 

006620 

0.007406 

2 

0 

002673 

0.004274 

0,005313 

0, 

006282 

0 

008379 

0,009240 

3 

0 

003671 

0,005508 

0*006668 

0, 

007735 

0 

010010 

0,010934 

4 

0 

004669 

0.006712 

0,007978 

0, 

009130 

0 

011561 

0,012540 

5 

0 

005668 

0.007894 

0,009255 

0, 

010484 

0 

013055 

0,014085 

6 

0 

006667 

O.OO 9 O 6 I 

0,010508 

0. 

011808 

0 

014508 

0,015584 

, 7 

0 

007667 

0,010216 

0,011743 

0, 

013108 

0 

015928 

0,017047 

a 

0 

008666 

0.011361 

0,012962 

0, 

014388 

0 

017321 

0,018480 

9 

0 

009665 

0,012497 

0,014169 

0. 

015653 

0 

018691 

0,019889 

10 

0 

010665 

0,013626 

0,015365 

0. 

016903 

0 

020043 

0,021276 

11 

0 

011665 

0,014749 

0,016552 

0, 

018142 

0 

021377 

0,022645 

12 

0 

012664 

0.015866 

0.017730 

0. 

019370 

0 

022697 

0,023999 

13 

0 

013664 

0,016979 

0,018901 

0, 

020589 

0 

024005 

0.025338 

14 

0 

014663 

0,018088 

0,020066 

0* 

021800 

0 

025301 

0*026664 

15 

0 

015663 

0.019192 

0,021225 

0, 

023004 

0 

026586 

0,027979 

16 

0 

016662 

0.020294 

0»022379 

0* 

024200 

0 

027862 

0,029283 

17 

0 

017662 

0.021392 

0.023528 

0, 

025391 

0 

029129 

0,030578 

18 

0 

018662 

0,022487 

0.024673 

0, 

026576 

0 

030389 

0,031865 

19 

0 

019661 

0,023579 

0,025813 

0, 

027755 

0 

031641 

0,033143 

20 

0 

020661 

0,024669 

0,026950 

0, 

028930 

0 

032886 

0,034414 

21 

0 

021660 

0,025757 

0,028083 

0, 

030101 

0 

034125 

0,035678 

22 

0 

022660 

0,026843 

0,029213 

0. 

031267 

0 

035358 

0.036935 

23 

0 

023660 

0,027926 

0,030340 

0, 

032429 

0 

036586 

0,038186 

24 

0 

024659 

0.029008 

0,031463 

0* 

033587 

0 

037809 

0.039432 

25 

0 

025659 

0,030088 

0,032585 

0, 

034743 

0 

039026 

0,040672 

26 

0 

026659 

0,031166 

0,033703 

0, 

035894 

0 

040239 

0,041907 

27 

0 

027658 

0,032242 

0,034820 

0, 

037043 

0 

041448 

0,043138 

26 

0 

028658 

0,033317 

0,035934 

0, 

038189 

0 

042653 

0,044364 

29 

0 

029658 

0,034391 

0,037045 

0. 

039332 

0 

043854 

0,045585 

30 

0 

030657 

0,035463 

0,038155 

0, 

040472 

0 

045050 

0,046803 

31 

0 

031657 

0,036534 

0,039263 

0, 

041610 

0 

046244 

0,048016 

32 

0 

032657 

0,037604 

0,040369 

0. 

042746 

0 

047434 

0,049226 

33 

0 

033656 

0,038673 

0.041473 

0, 

043879 

0 

048621 

0,050432 

34 

0 

034656 

0,039740 

0,042575 

0, 

045010 

0 

049804 

0,051635 

35 

0 

035655 

0,040806 

0.043676 

0, 

046138 

0 

050985 

0,052834 

36 

0 

036655 

0,041872 

0,044775 

0* 

047265 

0 

052163 

0,054030 

37 

0 

037655 

0,042936 

0,045873 

0, 

048390 

0 

053338 

0,055224 

38 

0 

038654 

0,043999 

0,046969 

0, 

049513 

0 

054511 

0.056414 

39 

0 

039654 

0,045062 

0,048064 

0. 

050634 

0 

055681 

0.057602 

40 

0 

040654 

0,046123 

0.049157 

0. 

051753 

0 

056848 

0,058786 

41 

0 

041653 

0,047184 

0,050249 

0. 

052871 

0 

058013 

0,059969 

42 

0 

042653 

0,048244 

0,051340 

0, 

053987 

0 

059176 

0,061148 

43 

0 

043653 

0*049303 

0,052429 

0, 

055102 

0 

060337 

0,062326 

44 

0 

044652 

0,050361 

0,053518 

0* 

056215 

0 

061495 

0,063501 

45 

0 

045652 

0,051418 

0,054605 - 

0, 

057326 

0 

062652 

0,064673 

46 

0 

046652 

0,052475 

0,055691 

Op 

058436 

0 

063806 

0,065844 

47 

0 

047651 

0,053531 

0,056776 

0-. 

059545 

0 

064959 

0,067012 

48 

0 

048651 

0,054586 

0,057860 

0, 

060653 

0 

066109 

0,068178 

49 

0 

049651 

0,055641 

0,058943 

0. 

061759 

0 

067258 

0,069342 

50 

0 

050650 

0,056695 

0,060025 

0. 

062863 

0 

068405 

0,070504 





Figure 6 

MAIN Program Basic Plow C^rt 
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Figure 7 

Detailed Plow Chart of Capture Maneuver Logic 























Figure 7 (Continued) 

Detailed Flow Chart of Capture Maneuver Logic 











Figure 7 (Continued) 

Detailed Plow Chart of Capture Maneuver Logic 
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Figure 7 (Continued) 

Detailed Flow Chart of Capture Maneuver Logic 
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Figure 7 (Continued) 

Detailed Flow Chart of Capture Maneuver Logic 
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Figure 7 (Continued) 

Detailed Flow Chart of Capture Maneuver Logic 
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Figure 7 (Continued) 

Detailed Plow Chart of Capture Maneuver Logic 
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Figure 7 (Concluded) 

Detailed Flow Chart of Capture Maneuver Logic 
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Figure 9 

Sample Problem Input Data - Pitch Axis 
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4 

53000. 


500. 
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5 
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33 
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34 
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35 
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39 
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31 
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33 

0.1 



COARSE INTEGRATION STEP SIZE(SEC) 


33 

0.001 



MINIMUM INTEGRATION STEP SIZE (SEC) 

34 



Figure 9 (Concluded) 

Sample Problem Input Data - Pitch Axis 
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-0.416 THETA ERROR/ALPHA SENSITIVITY 1ST STAGE 39 
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0. 120. 4. 99. 8. 82. 47 
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10 MEAN VALUE OF THRUST MISALIGNMENT( DEG ) VS. TIME(SEC) 52 

0. 0.090 2. 0.094 5. 0.089 53 

10. 0.082 15. 0.080 54 

10 STANDARD DEVIATION OF THRUST MISALIGNMENT ( DEG ) VS. TIME(SEC) 55 

0. 0.060 1. 0.060 2. 0.034 56 

5. 0.025 15. 0.025 57 

20 NO. HISTOGRAM CELLS, MIN AND MAX THETA ERROR, TIME, IMPULSE ( NHIST= 1 ) 58 

0. 10. 0. 10. 0. 2000. 59 


*EOR 



Figure 10 

Sample Problem Input Data - Yaw Axis 


I 


I 


a 0 199 

50 

3 3 

1 0 


1 

SAP1PLE PROBLEM 

SCOUT 

SECOND STAGE VAU CAPTURE NO INTEGRAL GAIN 


3 

199 SAMPLES, NPRINT-0, 

INC-50, 

3 OUERSHOOTS, 13 SEC. TIME LIMIT 


3 

1334567 



INITIAL RANDOM NUMBER SEED (K) 


4 

53000. 


500. 

MOMENT OF INERTIA (SLUG-FT3) 

MEAN, SIG 

5 

378.8 


1 . 

CENTER OF GRAUITY STATION (IN.) 

MEAN, SIG 

6 

1. 


0.0158 

BOOSTER THRUST RATIO TO NOM. 

MEAN, SIG 

7 

1. 


0.13 

Q RATIO TO NOMINAL 

MEAN, SIG 

8 

0.803 


0.0367 

DEADBAND HALFUIDTH (DEG) 

MEAN, SIG 

9 

0.035 


0.0167 

HYSTERESIS RATIO 

MEAN, SIG 

10 

0.073 


0.0088 

CONTROL MOTOR TURN-ON DELAY (SEC) 

MEAN, SIG 

11 

0 . 035 


0.0053 

CONTROL MOTOR TURN-OFF DELAY(SEC) 

MEAN, SIG 

13 

0.50 


0.0167 

RATE GAIN RATIO KR/KD(SEC) 

MEAN, SIG 

13 

494.1 


17.3 

RIGHT CONTROL MOTOR THRUST(LBS) 

MEAN, SIG 

14 

494.1 


17.3 

LEFT CONTROL MOTOR THRUST (LBS) 

MEAN, SIG 

15 

0.046 


0.0034 

CONTROL ELECTRONICS DELAY(SEC) 

MEAN, SIG 

16 

0.650 


0.0667 

RATE GYRO DAMPING FACTOR 

MEAN, SIG 

17 

33.5 


1 . 17 

RATE GYRO FREQUENCY (HZ) 

MEAN, SIG 

18 

-3.38 


0.897 

LOG OF ATT. OSCILLATION AMP. (DEG) 

MEAN, SIG 

19 

0.085 


0.196 

INITIAL RATE (DEG/SEC) 

MEAN, SIG 

30 

135. 


134. 

UIND UEL. AT IGNITION (FT/SEC) 

MEAN, SIG 

31 

148. 


153. 

UIND UEL. AT 3ND TIME (FT/SEC) 

MEAN, SIG 

33 

90. 


0.73 

YAU PATH ANGLE AT IGNITION ( DEG )MEAN, SIG 

33 

0.4735 


0.0037 

CNAS(FT3/DEG), A3(FT3/DEG3) 


34 

95. 


6.6 

XCP0(IN.), DXCP/DA( IN/DEG ) 


35 

467.68 



CONTROL MOTOR STATION (IN) 


36 

17. 



CONTROL MOTOR RADIAL LOCATION ( IN . ) 


37 

5. 



CONTROL MOTOR CANT ANGLE (DEG) 


38 

448.51 



BOOSTER NOZZLE THROAT LOCATION (IN 

. ) 

39 

1000.0 



BENDING FREQ. (HZ) (LARGE UALUE - SMALL DT ) 

30 

0. 



INITIAL TIME (SEC) 


31 

13. 



TIME LIMIT FOR CAPTURE(SEC) 


33 

0.1 



COARSE INTEGRATION STEP SIZE(SEC) 


33 

0.001 



MINIMUM INTEGRATION STEP SIZE (SEC) 

34 



Figure 10 (Concluded) 

Sample Problem Input Data - Yaw Axis 


I 

VJl 

0 

1 


8 


8 


6 

6 

20 


8 . 

4315. 

4995. 

0 . 

-0.416 

0 . 

10 . 

0 . 

NOMINAL DVNAMIC 

0 . 120 . 

12. 67. 


2ND TABLE TIME FOR UIND AND UELOCITV (SEC ) 
OELOCITY AT IGNITION (FT/SEC) 

UELOCITV AT 2ND TABLE TIME(FT/SEC) 
COMMANDED YAU PROGRAM RATE(DEG/SEC) 
THETA ERROR/ALPHA SENSITIUITY 1ST STAGE 
INITIAL ATTITUDE ERROR BIAS (DEG) 

ATT. ERROR TEST LIMIT FOR DIUERGENCE ( DEG ) 
NOMINAL ANGLE OF ATTACK BIAS (DEG) 
PRESSURECPSF ) US. TIME (SEC) 

4. 99. 8. 82. 


NOMINAL BOOSTER THRUST (LBS) US. TIME (SEC) 

0. 41000. 4.3 47155. 9.3 55936. 

12.5 61620. 

MEAN UALUE OF THRUST MISALIGNMENTC DEG ) US. TIME(SEC) 

0. 0.010 5. 0.003 10. 0.006 

STANDARD DEUIATION OF THRUST MISALIGNMENT ( DEG ) US. TIME(SEC) 

0. 0.048 5. 0.023 10. 0.020 

NO. HISTOGRAM CELLS, MIN AND MAX THETA ERROR, TIME, IMPULSE ( NHIST- 1 ) 

0. 10. 0. 10. 0. 2000. 


35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


*EOR 



Figure 11 

Sample Problem Input Bata - Boll Axis 


3 0 199 

50 

3 3 

1 0 


1 

SAMPLE PROBLEM 

SCOUT 

SECOND STAGE ROLL CAPTURE CASTOR I MOTOR 


3 

199 SAMPLES, NPRINT-3, 

INC-50, 

3 OUERSHOOTS, 13 SEC. TIME LIMIT 


3 

1334567 



INITIAL RANDOM NUMBER SEED (K) 


4 

435.1 


7.07 

MOMENT OF INERTIA (SLUG-FT3) 

MEAN.SIG 

5 

0 . 


0. 

:t**CARDS UITH ASTERISKS MUST BE AS SHOUN 

6 

1. 


0. 



7 

1 . 


0.13 

Q RATIO TO NOMINAL 

MEAN,SIG 

8 

1.433 


0.0476 

DEADBAND HALFUIDTH (DEG) 

MEAN,SIG 

9 

0.035 


0.0167 

HYSTERESIS RATIO 

MEAN,SIG 

10 

0.035 


0.0048 

CONTROL MOTOR TURN-ON DELAY (SEC) 

MEAN,SIG 

11 

0.016 


0.0030 

CONTROL MOTOR TURN-OFF DELAY(SEC) 

MEAN,SIG 

13 

0.45 


0.015 

RATE GAIN RATIO KR/KD(SEC) 

MEAN,SIG 

13 

90.90 


1.7 

RIGHT CONTROL MOTOR THRUST(LBS) 

MEAN,SIG 

14 

90.90 


1.7 

LEFT CONTROL MOTOR THRUST (LBS) 

MEAN,SIG 

15 

0.016 


0.0034 

CONTROL ELECTRONICS DELAY(SEC) 

MEAN,SIG 

16 

0.650 


0.0667 

RATE GYRO DAMPING FACTOR 

MEAN,SIG 

17 

36.5 


1.17 

RATE GYRO FREQUENCY (HZ) 

MEAN,SIG 

18 

-3.63 


0.66 

LOG OF ATT. OSCILLATION AMP. (DEG) 

MEAN,SIG 

19 

-0.04 


0.168 

INITIAL RATE (DEG/SEC) 

MEAN,SIG 

30 

0. 


0. 



31 

0. 


0. 



33 

0. 


0. 



33 

0. 


0. 

CLPHI*S(FT3/DEG ), A3 ( FT:t:|:3/DEG**3 

) 

34 

-31. 


0. 

-D REF (IN.) AERO REF. LENGTH 


35 

0. 





36 

17. 



ROLL MOTOR ARM (INCHES) 


37 

90. 





38 

687.6 





39 

100000. 



TORSIONAL FREQUENCY (HERTZ) 


30 

0. 



INITIAL TIME (SEC) 


31 

13. 



TIME LIMIT FOR CAPTURE (SEC) 


33 

0.1 



COARSE INTEGRATION STEP SIZE(SEC) 


33 

0.001 



MINIMUM INTEGRATION STEP SIZE (SEC) 

34 



Figure 11 (Concluded) 

Sample Problem Input Data - Roll Axis 


8 


4 

16 


16 


40 


1 . 

1 . 

1 . 

0 . 

0 . 

0 . 

10 . 

0 . 

NOniNAL DYNAHIC 

0. 1E0. 

la. 67. 


txt 

COPlflANDED ROLL PROGRAM RATE ( DEG/SEC ) 

*** 

INITIAL ATTITUDE ERROR BIAS (DEG) 

ATT. ERROR TEST LIMIT FOR DIOERGENCE ( DEG ) 
tXt 

PRESSURECPSF ) US. TIME (SEC) 

4. 99. 8. 8a. 


0 . 1 . 1000 . 

MEAN UALUE OF BOOSTER 

0 . 0 . 0 .a 


ttt TABLE NECESSARY FOR ROLL TORQUE 

1 . 

ROLL TORQUE (FT-LBS) US. TIME (SEC) 
-45. 0.4 -55. 


0.6 

-50. 

1.0 

4.0 

-10. 

15.0 

STANDARD DEUIATION OF 

0. 

a. 

0 .a 

0.6 

19. 

1.0 

4.0 

5. 

15.0 

NO. 

CELL FOR 

HISTOGRAM 

0. 

10. 

0. 


— a5 . 1.5 — 15. 

-5.0 

BOOSTER ROLL TORQUE ( FT-LBS ) US. TIME (SEC) 
a0. 0.4 E5. 

10. 1.5 8. 

4. 

MIN AND MAX THETA ERROR, TIME, IMPULSE 
10. 0. 3000. 


35 

36 

37 

38 

39 

40 

41 

4a 

43 

44 

45 

46 

47 

48 

49 

50 

51 

5a 

53 

54 

55 

56 

57 



Figure 12 

Group 1 - SPRINT = 0 


RUN NO. 1 PAGE NO. B 


SAMPLE PROBLEM SCOUT : 

SECOND STAGE 

PITCH UITH 

INTEGRAL 

GAIN 4B IN 

H/S 

199 SAMPLES, 

NPRINT-0, 

INC-50, 3 OUERSHOOTS, IB 

SEC TIME 

LIMIT, NHIST-1 

ITER 

MEAN 

SIGMA 

MEAN 

SIGMA 

MEAN 

SIGMA 

NO 

THETA 

THETA 

TIME 

TIME 

IMPULSE 

IMPULSE 

50 

1.9B199 

.B4179 

4.B6861 

1 .B144B 

756.96489 

B56. 15586 

100 

1.89088 

.B7689 

4. 144B3 

1 .39B58 

736.03906 

287.6B953 

150 

1.888B4 

.B58B8 

4.B44B5 

1 .B7771 

771.09885 

B81 .06BB9 

ERROR 

MODE 3ENC0UNTERED - 

DIUERGENT ATTITUDE UITH FOLLOUING 

data: 


1 

U1 

TIME (SEC)- 

3.9 

RATE(DEG/SE0- 

6.68 

ATT. ERROR (DEG) = 

10.309 

1 

Q 

(PSF) 

113. B 

ALPHA( DEGREE)- 

13. BB 

THRUST MIS(DEG)- 

.117 


F 

TOP (LB) - 

459. B 

F LOUER (LB) - 

470.81 

LAST OUERSHOOT - 

1.093 


199 

1.8977B 

.B5663 4.B9096 

1 .B0B81 

78B.31B30 B66. 

33504 



Figure 13 

Group 1 “ NPRINT “ 1 


RUN NO. 2 


PAGE NO. 2 


SAMPLE PROBLEM SCOUT 

SECOND STAGE 

VAU CAPTURE 

NO INTEGRAL GAIN 


199 SAMPLES, 

NPRINT-1, 

INC- 5, 3 OUERSHOOTS, 12 

SEC. TIME 

LIMIT 


ITER 

MEAN 

SIGMA 

MEAN 

SIGMA 

MEAN 

SIGMA 

NO 

THETA 

THETA 

TIME 

TIME 

IMPULSE 

IMPULSE 

5 

.72590 

.11674 

1.38240 

.77183 

147.99451 

164.80215 


THETA 

THETADOT 

TIME 

IMPULSE 


(DEG ) 

(DEG/SEC ) 

(SEC) 

(LB-SEC )/ 


.5445 

0.0000 

.4457 

0 . 0000 


-.7093 

0.0000 

2.6655 

145.2359 

(Case 6 - Overshoot Data 

-.6621 

0.0000 

2.9367 

170.5716 



THETA 

THETADOT 

TIME 

IMPULSE 


(DEG) 

(DEG/SEC ) 

(SEC) 

(LB-SEC )/ 


.4733 

0.0000 

1.1606 

95.1301 


.4673 

0.0000 

1.2619 

103.3494 

(Case 7) 

.6446 

0 . 0000 

1.8062 

149.1225 


THETA 

THETADOT 

TIME 

IMPULSE 


(DEG) 

(DEG/SEC ) 

(SEC) 

(LB-SEC )/ 


.6849 

0 . 0000 

1.4954 

302.9446 

(Case 8) 

.6445 

0 . 0000 

1.7040 

334.9332 

.7504 

0 . 0000 

2.0416 

386.7701 



THETA 

THETADOT 

TIME 

IMPULSE 

(DEG) 

(DEG/SEC ) 

(SEC) 

(LB-SEC ) 



Figure 14 

Group 1 - NPRINT = 2 


RUN NO. 

2 




PAGE 

NO. 2 

SAMPLE PROBLEM SCOUT 
E5 SAMPLES, NPRINT-E, 

SECOND STAGE VAU CAPTURE 
INC- 5, 3 OUERSHOOTS, 12 

NO INTEGRAL GAIN 
SEC. TIME LIMIT 


ITER 

NO 

MEAN 

THETA 

SIGMA 

THETA 

MEAN 

TIME 

SIGMA 

TIME 

MEAN 

IMPULSE 

SIGMA 

IMPULSE 

5 

.72590 

.11674 

1.38240 

.77183 

147.99451 

164.80215 


THETA 

THETADOT 

TIME 

IMPULSE 


(DEG ) 

(DEG/SEC ) 

(SEC) 

(LB-SEC )/ 


.5445 

0 . 0000 

.4457 

0.0000 


-.2402 

-1.2546 

1.8594 

0.0000 

(Case 6 - Critical Points 

-.4014 

-1.4209 

1.9799 

0.0000 

Deadband Crossing, 

-.5239 

-.6164 

2.1002 

57.4454 

Thrust -on & off, 

-.5524 

-.0355 

2.1874 

99.0876 

Overshoot points . ) 

-.6140 

-.4400 

2.4483 

99.0876 

1 

-.6787 

-.6340 

2.5688 

99.0876 


-.7055 

-.2241 

2.6313 

128.9294 


-.7093 

0 . 0000 

2.6655 

145.2359 


-.7001 

.3482 

2.7186 

170.5716 


-.6621 

0.0000 

2.9367 

170.5716 



THETA 

THETADOT 

TIME 

IMPULSE 

(DEG ) 

(DEG/SEC ) 

(SEC) 

(LB-SEC)/ 

.1847 

1.2144 

.8531 

0.0000 

.3485 

1,3782 

.9795 

0 . 0000 

.4529 

.5576 

1.0873 

56.6175 (Case 7) 

.4733 

0.0000 

1.1606 

95.1301 

.4724 

-.1190 

1 . 1763 

103.3494 

.4673 

0.0000 

1.2619 

103.3494 



Figure 15 

Group 1 - NPRINT = 3 


SAMPLE PROBLEM 

SCOUT 

SECOND STAGE 

VAU CAPTURE 

NO INTEGRAL GAIN 


as SAMPLES, NPRINT-3, 

INC« 5, 3 OUERSHOOTS, 13 

SEC. TIME 

LIMIT 


ITER 

MEAN 

SIGMA 

MEAN 

SIGMA 

MEAN 

SIGMA 

NO 

THETA 

THETA 

TIME 

TIME 

IMPULSE 

IMPULSE 

10 

.71774 

.08504 

1.67518 

.66616 

195.37197 

156.8E787 


THETA 

THETADOT 

TIME 

IMPULSE 

(DEG ) 

(DEG/SEC ) 

(SEC ) 

(LB-SEC )/ 

.4938 

0 . 0000 

.0307 

0.0000 

. 4915 

-.0380 

. 1000 

0.0000 (Case 11 - all integration 

.4849 

-.0934 

.3000 

0.0000 steps) 

. 4738 

-.1497 

.3000 

0 . 0000 

.4550 

-.3073 

.4000 

0.0000 

. 4313 

-.3664 

.5000 

0 . 0000 

.4016 

-.3379 

.6000 

0.0000 

.3656 

-.3933 

.7000 

0 . 0000 

.3330 

-.4597 

.8000 

0 . 0000 

.3734 

-.5311 

. 9000 

0.0000 

.3165 

-.6068 

1 . 0000 

0 . 0000 

. 1518 

-.6875 

1 . 1000 

0.0000 

.0787 

-.7737 

1.3000 

0 . 0000 

-.0033 

-.8663 

1.3000 

0 . 0000 

-.0948 

-.9656 

1.4000 

0 . 0000 

-.1968 

-1.0736 

1.5000 

0 . 0000 

-.3339 

-1 . 1006 

1.5350 

0 . 0000 

-.3378 

-1.1147 

1.5375 

0 . 0000 

-.3448 

-1 . 1319 

1.5438 

0 . 0000 

-.3483 

-1.1355 

1.5469 

0 . 0000 

-.3491 

-1 . 1364 

1.5477 

0.0000 

-.3917 

-1.3701 

1.6667 

0.0000 



Figure 15 (Continued) 
Group 1 - NPfilNT = 5 


I 


-.4464 

-.9146 

1.7167 

24.1102 

-.4670 

-.7376 

1.7417 

36.1653 

-.4757 

-.6493 

1.7542 

42.1928 

-.4796 

-.6051 

1.7604 

45 . 2066 

-.4805 

-.5941 

1.7620 

45.9600 

-.4810 

-.5886 

1.7628 

46.3367 

-.4815 

-.5831 

1 .7635 

46.7135 

-.5055 

0.0000 

1.8461 

86.5209 

-.5053 

.0538 

1.8537 

90.1907 

-.5042 

0.0000 

1.8947 

90.1907 


THETA 

THETADOT 

TlflE 

IfIPULSE 

(DEG ) 

(DEG/'SEC ) 

(SEC ) 

(LB-SEC )/ 

-2.2569 

.8738 

. 1284 

0 . 0000 

-2.1106 

2.0518 

.2284 

48.4547 

-2.0556 

2.3490 

.2534 

60.5683 

-2.0407 

2.4235 

.2596 

63.5967 

-2.0330 

2.4607 

.2628 

65.1110 

-2.0311 

2.4701 

.2635 

65.4895 

-1.7096 

3.7277 

.3673 

115.7545 

-1.3169 

4.1279 

.4673 

115.7545 

-1.2649 

4.1801 

.4798 

115.7545 

-1.2617 

4.1833 

.4806 

115.7545 

-.6874 

4.7622 

.6089 

115.7545 

-.2298 

4.3906 

.7089 

165.8308 

. 1928 

4.0610 

.8089 

215.9071 

.5846 

3.7755 

.9089 

265.9835 

.9500 

3.5317 

1.0089 

316.0598 

1.2929 

3.3276 

1 . 1089 

366.1361 

1.6174 

3.1615 

1.2089 

416.2124 

1.9271 

3.0318 

1.3089 

466.2887 

2.2255 

2.9374 

1.4089 

516.3650 

2.5163 

2.8776 

1.5089 

566.4413 

2.8027 

2.8520 

1.6089 

616.5176 


(Case 12) 



Figure 15 (Concluded) 




Group 1 - NPRINT 

■= 3 

3»0884 

2.8606 

1.7089 

666.5939 

3.3766 

2.9038 

1.8089 

716.6702 

3.6709 

2.9827 

1.9089 

766.7465 

3.9750 

3.0984 

2.0089 

816.8229 

4 .8935 

3.2530 

2.1089 

866.8992 

4.6276 

3.4488 

2.2089 

916.9755 

4.9845 

3.6888 

2.3089 

967.0518 

5.3678 

3.9768 

2.4089 

1017.1281 

5.7825 

4.3170 

2.5089 

1067.2044 

6.2340 

4.7143 

2.6089 

1117.2807 

6.7285 

5.1745 

2.7089 

1167.3570 

7.2724 

5.7037 

2.8089 

1217.4333 (Case 12 - Divergent) 

7.8730 

6.3086 

2.9089 

1267.5096 

8.5382 

6.9960 

3.0089 

1317.5860 

9.2766 

7.7722 

3.1089 

1367.6623 

ERROR nODE 3 

ENCOUNTERED - DIUERGENT 

ATTITUDE UITH FOLLOUING DATA: 

TIME (SEC) 

> 3.1 

RATE(DEG/SEO- 

7.77 ATT. ERROR(DEG)- 10 

Q (PSD 

« 112.6 

ALPHA( DEGREE)- 

15.09 THRUST niS(DEG)- 


F TOP(LB) 

= 500.8 

F LOUER (LB) 

- 484.55 LAST OUERSHOOT 

THETA 

THETADOT 

TIME 

IMPULSE 

(DEG ) 

(DEG/SEC ) 

(SEC) 

(LB-SEC )/ 

-1.2350 

.2441 

.1305 

0 . 0000 

-1.2127 

.6508 

.1805 

24.5905 

-1.2039 

.7526 

.1930 

30.7381 

-1 . 1990 

.8035 

.1992 

33.8120 , , 

-1 . 1978 

.8163 

.2008 

34.5804 (Case 13) 

-1 . 1971 

.8226 

.2016 

34.9646 

-1.0917 

1.5524 

.2903 

78.6147 

-.9378 

1.5257 

.3903 

78.6147 


097 

011 

000 



♦ 


Figure 16 

Group 2 Output - Unarranged Results 


RUN NO. 2 


SANPLE PROBLEN SCOUT 

SECOND STAGE 

YAU CAPTURE 

25 SANPLES, NPRINT-0, 

INC- 5, 3 OUERSHOOTS, 12 

CAP. ANGLE 

CAP. TINE 

INPULSE 

(DEGREES) 

(SECONDS) 

(LB-SEC) 

.85282 

1.32139 

406.56157 

.750G9 

2.15986 

11.77130 

.70802 

1.35056 

209.85397 

.53957 

1.91740 

90.77326 

.77838 

.16278 

21.01247 

.70933 

2.66549 

145.23594 

.64459 

1.80620 

149.12251 

.75038 

2.04156 

386.77013 

.76496 

1.50510 

424.64849 

.67867 

1 .82151 

107.97003 

.50552 

1.84609 

86.52092 

.74647 

2.94102 

297.10001 

.67147 

2.28236 

317.82994 

.87586 

.01637 

0.00000 

1.43315 

.06736 

0.00000 

.71359 

1.94375 

157.96780 

.75209 

2.40397 

289.30930 

.65715 

1.89176 

41.18371 

.73182 

2.91294 

178.54109 

1.30457 

.02874 

0.00000 

.68130 

1.81563 

126.37307 

.98769 

.01882 

0.00000 

.70236 

3.21081 

113.31923 

.70297 

3.98907 

161.24542 

.78631 

2.11634 

306.76177 


- NHIST = 0 


PAGE NO. 3 

NO INTEGRAL GAIN 
SEC. TINE LINIT 



Figure 17 

Group 3 Output - Histogram Data - NHIST •= 1 


RUN NO. 1 PAGE NO. 

SAMPLE PROBLEM SCOUT SECOND STAGE PITCH UITH INTEGRAL GAIN 42 IN H^S 
199 SAMPLES. NPRINT-O, INC-50, 3 OUERSHOOTS, 12 SEC TIME LIMIT, NHIST- 1 


CAPTURE ANGLE NO TIME NO IMPULSE NO 


0.00 

TO 

.50, 

0 

0 . 00 

TO 

.50, 

15 

0.00 

TO 

100.00, 

15 

.50 

TO 

1.00, 

2 

.50 

TO 

1.00, 

0 

100 . 00 

TO 

200.00, 

0 

1.00 

TO 

1.50, 

15 

1.00 

TO 

1.50, 

0 

200 . 00 

TO 

300.00, 

1 

1.50 

TO 

2.00, 

108 

1.50 

TO 

2.00, 

1 

300 . 00 

TO 

400.00, 

1 

2.00 

TO 

2.50, 

73 

2.00 

TO 

2.50, 

0 

400 . 00 

TO 

500 . 00 , 

2 

2.50 

TO 

3.00, 

1 

2.50 

TO 

3.00, 

1 

500.00 

TO 

600.00, 

7 

3.00 

TO 

3.50, 

0 

3.00 

TO 

3.50, 

0 

600.00 

TO 

700.00, 

19 

3.50 

TO 

4.00, 

0 

3.50 

TO 

4.00, 

0 

700.00 

TO 

800.00, 

45 

4.00 

TO 

4.50, 

0 

4.00 

TO 

4.50, 

43 

800.00 

TO 

900.00, 

45 

4.50 

TO 

5.00, 

0 

4.50 

TO 

5.00, 

136 

900.00 

TO 

1000.00, 

32 

5.00 

TO 

5.50, 

0 

5.00 

TO 

5.50, 

3 

1 000 . 00 

TO 

1100.00, 

21 

5.50 

TO 

6 . 00 , 

0 

5.50 

TO 

6.00, 

0 

1100.00 

TO 

1200.00, 

10 

6.00 

TO 

6.50, 

0 

6.00 

TO 

6.50, 

0 

1200.00 

TO 

1300.00, 

0 

6.50 

TO 

7.00, 

0 

6.50 

TO 

7.00, 

0 

1300.00 

TO 

1400.00, 

1 

7.00 

TO 

7.50, 

0 

7.00 

TO 

7.50, 

0 

1400.00 

TO 

1500.00, 

0 

7.50 

TO 

8.00, 

0 

7.50 

TO 

8.00, 

0 

1500.00 

TO 

1600.00, 

0 

8.00 

TO 

8.50, 

0 

8.00 

TO 

8.50, 

0 

1600.00 

TO 

1700.00, 

0 

8.50 

TO 

9.00, 

0 

8.50 

TO 

9.00, 

0 

1700.00 

TO 

1800.00, 

0 

9.00 

TO 

9,50, 

0 

9.00 

TO 

9.50, 

0 

1800.00 

TO 

1900.00, 

0 

9.50 

TO 

10.00, 

0 

9.50 

TO 

10.00, 

0 

1 900 . 00 

TO 

2000.00, 

0 



Figure 18 

Group 4 Output - l^iscrete Probability 


RUN NO. 1 PAGE NO. 4 

SAMPLE PROBLEM SCOUT SECOND STAGE PITCH UITH INTEGRAL GAIN 43 IN H/S 
199 SAMPLES, NPRINT-0, INC-50, 3 OUERSHOOTS, 13 SEC TIME LIMIT, NHIST- 1 


PROB CAP. ANGLE CAP. TIME IMPULSE 

(DEGREES) (SECONDS) (LB-SEC) 


.005 

.95337 

. 16011 

6.13147 

.010 

.99407 

. 16515 

31.19936 

.015 

1.03178 

.17519 

31.31533 

.030 

1.17653 

. 18857 

33.60156 

.035 

1 . 19316 

. 19801 

30.90811 

.030 

1.30856 

.30157 

33.40953 

.035 

1.31071 

.31497 

40.31509 

.040 

1.38335 

.31898 

43.33989 

.045 

1.38845 

.33300 

46.15758 

.050 

1.39451 

.33487 

46.73905 

.055 

1.39831 

.33870 

50.65938 

.060 

1.36888 

.33016 

51.53070 

.065 

1.37766 

.33849 

53.33380 

.070 

1.38374 

.36741 

57.59848 

.075 

1.44333 

.41766 

71.45598 

.080 

1.47849 

1.63489 

387.30133 

.085 

1 .49196 

3.83013 

331.79393 

.090 

1.50558 

4.35585 

433.70693 

.095 

1.51309 

4.31106 

460.44663 

. 100 

1.54064 

4.31537 

501.71154 

. 105 

1.54156 

4.34534 

534.59690 

.110 

1.56333 

4.35449 

538.11011 

.115 

1.57396 

4.36009 

543.53905 

. 130 

1.60986 

4.37473 

544.36847 

. 135 

1 .61019 

4.37773 

559.79617 



Figure 18 (Continued) 

Group 4 Output - Discrete Probability 


RUN NO. 1 PAGE NO. 

SAMPLE PROBLEM SCOUT SECOND STAGE PITCH UITH INTEGRAL GAIN 4E IN H/S 
199 SAMPLES, NPRINT-0, INC-50, 3 OUERSHOOTS, 12 SEC TIME LIMIT, NHIST- 1 


PROB 

CAP. ANGLE 
(DEGREES ) 

CAP. TIME 
(SECONDS) 

IMPULSE 
(LB-SEC ) 

.240 

1.81497 

4.46498 

706.62076 

.245 

1.81838 

4.47001 

707.17155 

.250 

1.82636 

4.47665 

707.97991 

.255 

1.82671 

4.48028 

709.46333 

.260 

1 .82762 

4.48157 

712.70609 

.265 

1.83541 

4.48232 

714.90028 

.270 

1.83822 

4.48477 

717.69771 

.275 

1.84410 

4.48700 

720.03541 

.280 

1.84778 

4.48779 

720.04334 

.285 

1.86167 

4.49045 

720.36123 

.290 

1.86681 

4.49172 

722.20390 

.295 

1.86695 

4.49416 

732.08493 

.300 

1 .87215 

4.49798 

732.33217 

.305 

1.87320 

4.50073 

738.46872 

.310 

1.87714 

4.50289 

739.32199 

.315 

1.87875 

4.50618 

740.14204 

.320 

1.87891 

4.51013 

743.81406 

.325 

1 .88198 

4.51105 

748.82072 

.330 

1.88405 

4.51246 

749.38205 

.335 

1.88504 

4.51436 

753.44694 

.340 

1.88680 

4.51582 

753.62911 

.345 

1.88733 

4.51816 

753.68218 

. 350 

1 .89251 

4.52320 

755.55886 

.355 

1.90160 

4.52447 

761.80032 

.360 

1.90191 

4.52605 

762.06797 



Figure 18 (Continued) 

Group 4 Output - Blscrete Probability 


RUN NO. 1 PAGE NO. 

SAMPLE PROBLEM SCOUT SECOND STAGE PITCH UITH INTEGRAL GAIN 43 IN H/S 
199 SAMPLES, NPRINT-0, INC-50, 3 OUERSHOOTS, 13 SEC TIME LIMIT, NHIST- 1 


PROB 

CAP. ANGLE 
(DEGREES ) 

CAP. TIME 
(SECONDS) 

IMPULSE 
(LB-SEC ) 

.475 

1.95174 

4.58403 

810.93177 

.480 

1.95331 

4.58636 

813.69659 

.485 

1.95390 

4.59933 

813.58913 

.490 

1 .95531 

4.60059 

831.43338 

.495 

1.95575 

4.60377 

833.83933 

.500 

1.95631 

4.60384 

835.57773 

.505 

1.96469 

4.60431 

837.10577 

.510 

1.96863 

4.60683 

837.83367 

.515 

1.96939 

4.60745 

831 . 13156 

.530 

1.96993 

4.61147 

833.88485 

.535 

1 .96994 

4.61934 

834.68096 

.530 

1.97147 

4.61989 

834.83876 

.535 

1.97344 

4.63936 

835.63935 

.540 

1.97383 

4.63885 

835.70038 

.545 

1.97389 

4.64553 

836.04854 

.550 

1.97433 

4.64800 

839.10901 

.555 

1.97658 

4.65331 

849.78964 

.560 

1.97784 

4.66438 

854.47960 

. 565 

1.97985 

4.66806 

857.70405 

.570 

1 .98043 

4.68438 

860.35033 

.575 

1 .98166 

4.68637 

860.61986 

.580 

1.98316 

4.68754 

861 .84803 

.585 

1.98533 

4.68755 

863.57430 

.590 

1.98836 

4.68908 

866.05315 

.595 

1.98855 

4.68933 

868.64973 



Figure 18 (Continued) 

Group 4 Output - Discrete Probability 


RUN NO. 1 PAGE NO. 

SAMPLE PROBLEM SCOUT SECOND STAGE PITCH UITH INTEGRAL GAIN 43 IN H/S 
199 SAMPLES, NPRINT-0, INC-50, 3 OUERSHOOTS, 13 SEC TIME LIMIT, NHIST- 1 


PROB 

CAP. ANGLE 
(DEGREES ) 

CAP. TIME 
(SECONDS) 

IMPULSE 

(LB-SEC) 

.710 

3.03394 

4.74556 

913.58711 

.715 

3.03445 

4.74678 

933.74508 

.730 

3.03485 

4.74798 

937.31875 

.735 

3.03714 

4.75103 

937.89049 

.730 

3.03778 

4.75399 

940.05338 

.735 

3.03796 

4.75694 

940.67391 

.740 

3.03933 

4.75706 

943.93587 

.745 

3.04340 

4.75896 

943.66154 

.750 

3.04885 

4.76133 

943.94383 

.755 

3.04980 

4.76178 

944.01645 

.760 

3.05369 

4.76313 

945.61506 

.765 

3.05364 

4.77636 

948.36703 

.770 

3.05498 

4.77960 

950.55707 

.775 

3.05756 

4.78331 

951.84339 

.780 

3.05957 

4.78350 

957.14790 

.785 

3.05993 

4.78373 

957.17877 

.790 

3.06014 

4.78783 

958.43993 

.795 

3.06079 

4.78855 

961,63735 

.800 

3.06933 

4.78918 

963.07453 

.805 

3.07095 

4.79054 

965.45803 

.810 

3.07113 

4.79587 

978.05693 

.815 

3.07318 

4.80019 

980,30018 

.830 

3.07334 

4.80135 

984.11719 

.835 

3.07733 

4.80130 

985. 17336 

.830 

3.08161 

4.80631 

986.93587 

.835 

3.08195 

4.80734 

997.68371 



Figure 18 (Concluded) 

Group 4 Output - Id.screte Probability 


SAMPLE PROBLEM SCOUT SECOND STAGE PITCH UITH INTEGRAL GAIN 43 IN H/S 
199 SAMPLES, NPRINT-0, INC-50, 3 OUERSHOOTS, 13 SEC TIME LIMIT, NHIST-1 


PROS 

CAP. ANGLE 
(DEGREES) 

CAP. TIME 
(SECONDS) 

IMPULSE 

(LB-SEC) 

.945 

3.13970 

4.93310 

1106.88744 

.950 

3.14761 

4.93313 

1109.64838 

.955 

3.15164 

4.94334 

1114.15347 

.960 

3.17061 

4.94561 

1114.80366 

.965 

3.18335 

4.96196 

1130.37638 

.970 

3.18749 

4.96437 

1136.33876 

.975 

3.18941 

4.97130 

1180.31498 

.980 

3.33344 

4.97389 

1183.40544 

.985 

3.37390 

5.00607 

1183.19656 

.990 

3.33136 

5.11458 

1196.01366 

.995 

3.99837 

5.19169 

1338.36315 


o> 

vn 


I 



Figure 19 

Group 5 Output - Abnormal Termination Summary 


RUN NO. 1 PAGE NO 

SAMPLE PROBLEM SCOUT SECOND STAGE PITCH UITH INTEGRAL GAIN 4E IN H^S 
199 SAMPLES. NPRINT-0. INC-50, 3 OUERSHOOTS, 12 SEC TIME LIMIT, NHIST- 1 


TOTAL NUMBER OF DIUERGENT CAPTURES - 1 


TOTAL NUMBER OF ABNORMAL CASES DELETED: 

DEADBAND OUERSHOOT - 0 

TIME LIMIT WITHOUT DIUERGENCE OR CAPTURE - 0 


FINAL UALUE OF RANDOM SEQUENCE INTEGER 


5585969 



APPENDIX 


FORTRAN PROGRAM LISTING 

A complete FORTRAN source program listing is presented in the following 
pages. It starts with the MAIN routine and is followed by the four ( 4 ) 
subroutines arranged in alphabetical order. There are a total of 627 cards in 
the MAIN routine. The total program including subroutines contains 723 cards. 


A-1 



»DECK nAIN 1 

PROGRAM MAINC INPUT, OUTPUT, TAPE5-INPUT,TAPES-0UTPUT) 2 

C 3 

C THIS IS A MONTE CARLO ANALYSIS OF THE STAGING CAPTURE MANUEUER 4 

C FOR A REACTION JET ATTITUDE CONTROL UITH A SIMPLE DEADBAND. IT 5 

C INCLUDES THE EFFECTS OF SEPARATION DISTURBANCES, NON-LINEAR 6 

C AERODYNAMIC MOMENTS, THRUST MISALIGNMENTS, AND OFF-NOMINAL 7 

C TRAJECTORY, UEHICLE, AND CONTROL SYSTEM CHARACTERISTICS. 8 

C 9 

INTEGER CAHIST(150),THIST(150),HISTIP(150) 10 

DIMENSION ETMEAN(50 ),ETSIG(50), TABQ ( 50 ) , TABT ( 50 ) , U ( 20 ) , XM ( 19 ) , 11 

1 XSD(19 ),THES( 1000 ),TTCAP(1000),TIMP(1000),MDl(3),nD3(3),MD4(3) 12 

COMMON RXD,Q,K,NPAGE,NLINE,NRUN,NTITLE(26) 13 

DATA MEl/l/,ME3/3/,ME4/4/,MDl(l)/6HDEADBA/,MDl(2)/6HND OUE/, 14 

1 MD1(3)/6HRSH00T/,MD3(1)/6HDIUERG/,MD3(2)/6HENT AT/,MD3(3) 15 

a /6HTITUDE/,MD4( 1 )/6HTIME L/, MD4 ( 2 )/6HIMIT B/, MD4 ( 3 )/6HENIGN / 16 

C READ IN INTEGER CONSTANTS AND RUN OPTIONS 17 

10 READC 5,1020) NRUN,NPRINT,NTRUN, INC,NSTAGE,N,NHIST,NKI 18 

C TEST INPUT FILE FOR END OF FILE, IF NO MORE DATA STOP. 19 

IF (EOF(5).NE.0) STOP 20 

C READ TUO CARDS OF HOLLERITH IDENTIFICATION OF RUN FOR OUTPUT 21 

READ( 5,1030) CNTITLE( J), J-1,26) 22 

IF(EOF(5).NE.0) STOP 23 

C INITIALIZE CONSTANTS AND SUMS 24 

STS-0. 25 

SSTS-0. 26 

STT“0. 27 

SSTT-0. 28 

STI-0. 29 

SSTI-0. 30. 

GINT-0. 31 

TKI-0. 32 

MQ-1 33 

MT-1 34 



O O o o 


f1En-l 35 

MES-l 36 

NPAGE-1 37 

NLINE-0 38 

NERRl-0 39 

NERR3-0 40 

NERR4-0 41 

C READ INITIAL RANDOM NUMBER GENERATOR SEED 42 

READC 5,1040) K 43 

C UARMUP PSEUDO-RANDOM NUMBER GENERATOR UITH 10 PASSES 44 

DO E0 J-1,10 45 

E0 CALL RNDX(l) 46 

C READ IN MEANS AND STANDARD DEUIATIONS OF RANDOM UARIABLES 47 

READ( 5,1050) (XM(I ),XSD(I ),!=■!, 19) 48 

C READ IN AERODYNAMIC COEFFICIENTS 49 

READ( 5,1050) CNAS, A3,XCP0,DCPDA 50 

C READ IN SINGLE UALUED CONSTANTS 51 

READ( 5,1060) XF,ZC,XLAM,XN,UB,TI,TMAX,DELT0,DELTMN,AT,UO,U1, 53 

1 GDOT, ASEN,THCMO,TLIMIT, ADES 53 

SL=SIN(XLAM/57.3 ) 54 

CL=C0S(XLAM/57.3 ) 55 

C TEST TO SEE IF INTEGRAL GAIN CONTROL IS USED. NKI-1 YES. 56 

IF (NKI .LE. 0 ) GO TO 30 57 

IF INTEGRAL GAIN USED, READ IN GAIN (KI), LIMITER UALUE XLIM, 58 

TIME AT UHICH INTEGRAL GAIN IS STARTED (TKI) 59 

READC 5,1060) GINT, XLIM, TKI 60 

READ IN TABLES OF DYNAMIC PRESSURE, THRUST, MEAN AND SIGMA THRUST 61 
MISALIGNMENT 63 

30 READ(5,1070) NTl , ( TABQ ( J ) , J • 1 , NTl ) 63 

READ<5,1070) NTS, (TABT(J),J-1,NT3) 64 

READ(5,1070) NT3, (ETMEAN( J ), J- 1,NT3 ) 65 

READ(5,1070) NT4 , ( ETSIG ( J ) , J - 1 , NT4 ) 66 

N3-INC 67 

NPH-0 68 
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J-0 69 

C TEST FOR HISTOGRAfI DATA, NHIST-1 IS VES 78 

IF (NHIST .EG. 0 ) GO TO 50 71 

C READ IN HISTOGRAM CELL INFORMATION 72 

READC 5,1070) NCEL,CMINH,CMAXH,TMINH,TMAXH,PMINH,PMAXH 73 

C SET UP HISTOGRAM CELL INFORMATION 74 

XNCEL-NCEL 75 

DO 40 LCEL-1,NCEL 76 

CAHIST(LCEL )-0 77 

THIST(LCEL)-0 78 

40 HISTIP(LCEL )-0 79 

DIFC»(CMAXH-CMINH)/XNCEL 80 

DIFT*(TMAXH-TMINH)/XNCEL 81 

DIFP-(PMAXH-PMINH)/XNCEL 82 

83 

THIS IS THE START OF THE MONTE CARLO LOOP, J IS THE SAMPLE NUMBER 84 

85 

50 J=J+1 86 

C COMPUTE UALUE FOR EACH OF THE NORMALLV DISTRIBUTED UARIABLES 87 

60 DO 70 1-1,16 88 

CALL RNDX(0) 89 

U(I )=XM(I) + RXD*XSD(I) 90 

70 CONTINUE 91 

RINT-U(l) 92 

CG=U(2) 93 

TDELTA=U(3) 94 

THEDD-0. 95 

QDELTA-U(4) 96 

DB-UC5) 97 

H-U(6) 98 

Tl-U(7) 99 

T2-U(8) 100 

GAIN-U(9) 101 

FT-U(10) 102 



FB-U(ll) 103 
FTAU-U(ia) 104 
GTAU-U( 13 14159*U( 14 ) ) 105 
THEDTO-U(IS) 106 
THED-0. 107 
CALL RNDX(l) 108 
ANGX-3. 14159*0-1 .57079 109 
CALL RNDXC0) 110 
GME=«RXD*XSD(19)-ADES 111 
CALL RNDX(0) 118 
RXDET=RXD 113 
COMPUTE RANDOM WIND AT TUO POINTS 114 
CALL RNDX(0) 115 
UU0=XM(17)+RXD*XSD(17) 116 
UU1=XM(18)+RXD*XSD(18) 117 
GM1-XM(19)+GD0T*AT 118 
COMPUTE ANGLE OF ATTACK DUE TO UIND AT IGNITION AND RATE OF CHANGE119 
AUPO-(57.3*UUO*SIN(XM(19)/57.3) )/(U0+UU0*C0S(XM(19)/57.3) ) 180 
AUPR-AUPO-GME 181 
AUPD= ( (57.3*UUl*SIN(GMl/57.3) )/(Ul+UUl*C0SCGMl/57.3) )-AUPO)/AT 188 
AUD-AUPD-GDOT 183 
CALL RNDX( 1 ) 124 
ANGV-3. 14159*0-1 .57079 185 
COMPUTE THE INITIAL ATTITUDE ERROR 126 
THE*EXP(U(15) )*SIN(ANGY )+AUPR*ASEN+THCMO 127 
BEGIN THE CAPTURE ANALYSIS 188 
SET ALL INITIAL CONDITIONS AND COMPUTE CONSTANTS 129 
NPK-0 130 
THETN-THE 131 
XINN-0. 132 
TIMPN-0.0 133 
TN-TI 134 
THEROR-THETN-GDOT*TN 135 
D-DB 136 
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RATI01-(XN-CG)/(ia.»RINT) 137 

THEDN-THED+THEDTO 138 

DBOFF-DB*( 1 .-H ) 139 

THECAP-0.0 140 

C COflPUT TURN-ON AND TURN-OFF TIME DELAYS 141 

TY-Tl+FTAU+GTAU 143 

TZ-Ta+FTAU+GTAU 143 

C COnPUTE CONTROL ACCELERATIONS FOR TOP AND BOTTOM MOTORS 144 

TDDCT = 4.775*FT*< (XF-CG )ScCL+ZC*SL )/RINT 145 

TDDCB--4.775*FB*( ( XF-CG )*CL+ZC*SL )/RINT 146 

C SET NO. OF OUERSHOOTS TO ZERO 147 

NOS=>0 148 

TEST NSTAGE. NSTAGE*a CONTROL SYSTEM IS ACTIOATED AT IGNITION 149 

NSTAGE=3 CONTROL SYSTEM IS ALREADY ACTIVATED 150 

IF (NSTAGE-a .LE. 0 ) GO TO 110 151 

C CONTROLS ALREADY ACTIUATED. FIRE MOTOR IF OUTSIDE DEADBAND 153 

THETS0=THEROR+GAIN*THEDN 153 

IF (ABS(THETS0)-D .LT. 0 ) GO TO 110 154 

C OUT OF DEADBAND 155 

IF (THETS0 .LT. 0 ) GO TO 80 156 

C OUT OF + DEADBAND, TURN ON BOTTOM MOTOR 157 

THEDDC=TDDCB 158 

FC=FB 159 

GO TO 90 160 

C OUT OF - DEADBAND, TURN ON TOP MOTOR 161 

80 THEDDC=TDDCT 163 

FC-FT 163 

C CHANGE DEADBAND TO TURN OFF LEUELCLESS HYSTERESIS) 164 

90 D-DBOFF 165 

GO TO 130 166 

C CONTROL ACTIUATED AT IGNITION, SET CONTROL FORCE TO ZERO 167 

110 THEDDC-0.0 168 . 

FC-0.0 169 

C COMPUTE EFFECTIUE DELAY TIME DUE TO BENDING, USE IF NON-NEGATIUE 170 
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1E0 TB-SIN(ANGX ) 171 

COMPUTE TOTAL TURN ON AND TURN OFF DELAY INCLUDING BENDING 173 

IF LESS THAN ZERO, SET DELAY TO ZERO 173 

TAUl-TB+TY 174 

IF (TAUl .GE. 0 ) GO TO 130 175 

TAUl-0. 176 

130 TAUa-TB+TZ 177 

IF (TAUa .GE. 0 ) GO TO 140 178 

TAUa-0. 179 

COMPUTE INITIAL ERROR SIGNAL 180 

140 THETS0=THEROR+GAIN*THEDN 181 

IF ( ABS(THETS0 )-D . LT , 0 ) GO TO 150 183 

OUT OF DEADBAND 183 

IF (THEDDC .NE. 0 ) GO TO 160 184 

OUT OF DEADBAND CONTROL MOTOR OFF 185 

GO TO 510 186 

150 IF (THEDDC . NE . 0 ) GO TO 510 187 

IN DEADBAND- CONTROL MOTOR OFF OR OUT OF DEADBAND CONTROL MOTOR 0N188 

SET INTEGRATION STEP SIZE BACK TO COARSE LEUEL 189 

160 DELTAT = DELT0 190 

C COMPUTE AUERAGE TIME FOR INTEGRATION AND COMPUTE ACCELERATIONS 191 

170 TN3=TN+DELTAT/a. 193 

CALL TBLU(NT1,0RD,TN3,TABQ,MQ ) 193 

DYNP»ORD*QDELTA 194 

CALL TBLU(NT3,0RD,TN3,TABT,MT) 195 

THRUST=ORD*TDELTA 196 

CALL TBLU ( NT3, ETM, TN3, ETMEAN, MEM ) 197 

CALL TBLU ( NT4, ETS, TN3, ETSIG, MES ) 198 

ET-ETM+RXDET*ETS 199 

C COMPUTE ANGLE OF ATTACK AND TOTAL ACCELERATION 300 

AL-AUPR+THETN+AUD^TN3+THEDN*DELTAT/3.+THEDD:«C0.5*DELTAT:«c»3 301 

TDDN3-THEDDC+RATI0lJt:ET:*:THRUST+4.775*DYNP*(CNAS*AL+A3*AL*AL*AL ) 303 

1»(CG-XCP0-DCPDA*ABS(AL ) )/RINT 303 

C SET TOTAL ACCELERATION 304 



A-8 


THEDD-TDDNa 305 

C COMPUTE RATE AT END OF INTEGRATION STEP 306 

THEDN1-TDDN3XDELTAT+THEDN 307 

C COMPUTE ATTITUDE AT END OF INTEGRATION STEP 308 

THETNl- (THEDNl+THEDN )*DELTAT/3 . +THETN 209 

C COMPUTE ATTITUDE ERROR AT END OF INTEGRATION STEP 210 

THER01-THETN1-GD0T^(TN+DELTAT ) 211 

C COMPUTE INTEGRAL OF ATTITUDE ERROR AT END OF INTEGRATION STEP 212 

XINNl=XINN+THEROR*DELTAT 213 

C COMPUTE INTEGRAL GAIN TERM AND LIMITED INTEGRAL GAIN FUNCTION 314 

EPIT1"GINT*XINN1 215 

EPIN1»SIGN(AMIN1(ABS(EPIT1),XLIM),EPIT1 ) 216 

C TEST TO SEE IF TIME TO START INTEGRAL GAIN LOGIC 217 

IF(TN.LE.TKI)EPIN1=0. 218 

C TEST FOR DIUERGENCE LIMIT ON ATTITUDE ERROR 219 

IF ( ABSCTHEROl )-TLIMIT .GE. 0 ) GO TO 440 220 

C STILL UITHIN ATTITUDE ERROR LIMITS- NO DIUERGENCE 221 

IF (THEDDC .EG. 0 ) GO TO 180 222 

C CONTROL MOTOR ON, COMPUTE ERROR SIGNAL AT END OF STEP 223 

THETS-THER01+GAIN*THEDN1+EPIN1 224 

IF (ABS(THETS)-D .LE. 0 ) GO TO 190 225 

C OUT OF DEADBAND 226 

GO TO 390 227 

180 THETS=THER01+GAIN*THEDN1+EPIN1 228 

IF (ABS(THETS)-D .LT. 0 ) GO TO 210 229 

C IN DEADBAND, TEST FOR MINIMUM INTEGRATION STEP SIZE 230 

190 IF (DELTAT-DELTMN .LE. 0 ) GO TO 210 231 

C CUT STEP SIZE IN HALF AND TRV AGAIN 232 

200 DELTAT=DELTAT/2. 233 

GO TO 170 234 

C TEST FOR ATTITUDE ERROR OUERSHOOT 235 

210 IF (THEDN-GDOT) 230,240,220 236 , 

220 IF (THEDNl-GDOT .LE. 0 ) GO TO 240 237 

C NO OUERSHOOT ENCOUNTERED ON THIS STEP 238 
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GO TO 330 S39 

E30 IF (THEDNl-GDOT .LT. 0 ) GO TO 330 240 

C POSITIUE RATE 241 

240 KCASE-0 242 

C COnPUTE OUERSHOOT ATTITUDE, RATE, ATTITUDE ERROR 243 

250 THET0S-(GD0T:*:*2-THEDN*)»c2)/(2.*TDDN2)+THETN 244 

THEDOS-GDOT 245 

THEEOS -THETOS-GDOTC ( TN+ ( GDOT-THEDN )/THEDD ) 246 

C UPDATE NO. OF OUERSHOOTS, IF INTEGRAL GAIN USED,ONLV AFTER START 247 

IFCTN.GE.TKDNOS-NOS+l 248 

C CONPUTE TINE AND INPULSE TO OUERSHOOT 249 

TX-( GDOT-THEDN )/TDDN2 250 

TOS-TN+TX 251 

TINPOS=TINPN+FC*TX 252 

C TEST OUERSHOOT UERSUS LAST OUERSHOOT FOR A NAXINUN 253 

IF (ABS(THEEOS)-THECAP .LT. 0 ) GO TO 260 254 

THECAP-ABS(THEEOS ) 255 

TINPCP-TINPOS 256 

TCAP=TOS 257 

260 IF(NPRINT.EQ.0) GO TO 310 258 

C TEST FOR PAGE CHANGE ON OUTPUT 259 

270 IF (NPK .NE. 0 ) GO TO 290 260 

280 URITE( 6,870) 261 

NLINE=NLINE+2 262 

URITE( 6,880) 263 

NLINE=NLINE+2 264 

NPK=NPK+1 265 

290 IF (NLINE+1-59 . LE . 0 ) GO TO 300 266 

NPAGE=-NPAGE + 1 267 

CALL PAGEHD 268 

GO TO 280 269 

C PRINT ATTITUDE ERROR, RATE, TINE, AND INPULSE AT OUERSHOOT 270 

300 URITE( 6,890) THEEOS, THEDOS, TOS, TINPOS 271 

NLINE-NLINE+1 272 
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C TEST TO SEE IF MAXinUM NUMBER OF OUERSHOOTS HAS BEEN REACHED E73 

310 IF (NOS-N .LT. 0 ) GO TO 3E0 274 

C THIS INDIUIDUAL CAPTURE RUN IS FINISHED 275 

GO TO 680 276 

C NOT ENOUGH OUERSHOOTS, CONTINUE 277 

320 IF (KCASE . GT . 0 ) GO TO 600 278 

C UPDATE IMPULSE AND TIME 379 

330 TIMPN-FC*DELTAT+TIMPN 280 

TN=TN+DELTAT 281 

IF (THEDDC .NE. 0 ) GO TO 340 282 

C MOTOR OFF, TEST FOR DEADBAND 283 

IF ( ABS(THETS )-D . GE , 0 ) GO TO 350 284 

GO TO 460 285 

340 IF (ABS(THETS )-D .GT. 0 ) GO TO 460 286 

C MOTOR OFF AND OUT OF DEADBAND, OR MOTOR ON AND IN DEADBAND 287 

350 THEDN'THEDNl 288 

THETN»THETN1 289 

THER0R*THER01 290 

XINN=XINN1 291 

C TEST PRINT OPTION 292 

IF (NPRINT-2 .LT. 0 ) GO TO 510 293 

IF (NPK .NE. 0 ) GO TO 370 294 

360 URITEC 6,870) 295 

URITEC 6,880) 296 

NLINE=NLINE+4 297 

NPK=NPK+1 298 

370 IF (NLINE+1-59 . LE . 0 ) GO TO 380 299 

NPAGE=NPAGE+1 300 

CALL PAGEHD 301 

GO TO 360 302 

C PRINT EACH POINT OF ATTITUDE ERROR, RATE, TIME, AND IMPULSE 303 

380 URITEC 6,890) THEROR, THEDN, TN, TIMPN 304 

GO TO 510 305 

390 IF (THEDDC .LT. 0 ) GO TO 400 306 
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C 


c 


c 


c 


IF (THETS .GE. 0 ) GO TO 410 307 

GO TO ai0 308 

400 IF (THETS .GT. 0 ) GO TO 310 309 

410 IF (DELTAT-DELTMN .GT. 0 ) GO TO 200 310 

IF (ABS(TN-TKI)-DELT0 .GT. 0 ) GO TO 430 311 

IF (ABS(EPINl) .EQ. 0 ) GO TO 430 31E 

IF (THETS .LT. 0 ) GO TO 420 313 

FC-FB 314 

THEDDC^TDDCB 315 

KCASE=0 316 

GO TO 160 317 

TOP MOTOR FIRING 318 

420 FC«FT 319 

THEDDC-TDDCT 320 

KCASE=0 321 

GO TO 160 322 

DETECT ERROR MODE 1 - DEADBAND OUERSHOOT 323 

430 NERRl-NERRl+1 324 


NLINE»NLINE+5 325 

URITE (6,1100) MEl , MDl ( 1 ) , MDl < 2 ) , MDl ( 3 ) , TN, THEDN, THEROl , DYNP, AL, ET, 326 
1 FB,FT,THECAP 327 

GO TO 60 328 


DIUERGENT TO LIMITING ATTITUDE ERROR , ERROR MODE 3 329 

440 NERR3=NERR3+1 330 

NLINE=NLINE+S 331 

URITE ( 6, 1 100 ) ME3, MD3 ( 1 ) , MD3 ( 2 ) , MD3 ( 3 ) , TN, THEDN, THEROl , DYNP, AL, ET, 332 
1 FB,FT,THECAP 333 

GO TO 60 334 

ERROR MODE 4 REACHED LIMITING TIME BEFORE ATTITUDE OUERSHOOT 335 

450 NERR4=NERR4+1 336 

NLINE-NLINE+5 337 

URITE(6, 1100) ME4,MD4(1 ), MD4 ( 2 ), MD4 C3 ), TN, THEDN, THEROl , DYNP, AL, ET, 338 


1 FB,FT,THECAP 


339 


IF (ABS( THEROl )-THECAP .LT. 0 ) GO TO 680 


340 
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C UPDATE CAPTURE DATA 341 

THECAP-ABS(THER01 ) 343 

TCAP-TN 343 

TIMPCP-TIMPN 344 

GO TO S80 345 

C UPDATE INTEGRATION UARIABLES 346 

460 THEDN=THEDN1 34? 

THETN'THETNl 348 

THEROR-THEROl 349 

XINN=XINN1 350 

C TEST PRINT OPTION 351 

IF (NPRINT-3 .NE, 0 ) GO TO 500 353 

IF (NPK .GT. 0 ) GO TO 480 353 

470 URITE( 6,870) 354 

URITEC 6,880) 355 

NLINE=NLINE+4 356 

NPK=NPK+1 357 

480 IF (NLINE+1-59 ,LE. 0 ) GO TO 490 358 

f NPAGE=NPAGE+1 359 

w CALL PAGEHD 360 

GO TO 470 361 

C PRINT ATTITUDE ERROR, RATE, TIME, AND TOTAL IMPULSE 363 

490 URITE( 6,890) THEROR, THEDN, TN, TIMPN 363 

C TEST FOR CASE TERMINATION DUE TO TIME LIMIT 364 

500 IF (TN-TMAX .LT. 0 ) GO TO 170 365 

REACHED TIME LIMIT. IF OUT OF DEADBAND CONTINUE UNTIL OUERSHOOT 0R366 
DIUERGENCE IS REACHED 367 

IF (ABS(THETS )-DB .GT. 0 ) GO TO 170 368 

C IN DEADBAND. IF ERROR IS NEAR LIMITING ERROR CONTINUE TO DIUERGENC369 

IF ( ABS(THEROR )-XLIM-DB .GT. 0 ) GO TO 170 370 

C IN DEADBAND. IF RATE IS LESS THAN MINIMUM MOTOR FIRING, TERMINATE 371 

IF (ABS(THEDN)-TDDCT^TAUS .LE. 0 ) GO TO 450 373 

C CONTINUE INTEGRATION 373 

GO TO 170 374 
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TEST FOR MOTOR ON OR OFF, AT SUITCHING LINE, SET INTEGRATION STEP 375 


TO ON OR OFF DELAY TIME. THIS LOGIC GOES TO STATEMENT 600 376 

510 IF (THEDDC ,EO. 0 ) GO TO 530 377 

C CONTROL MOTOR ON, SET INTEGRATION STEP TO TURN OFF DELAV(TAUa) 378 

DELTAT-TAUE 379 

GO TO 530 380 

C CONTROL MOTOR OFF, SET INTEGRATION STEP TO TURN ON DELAV(TAUl) 381 

530 DELTAT=TAU1 383 

C PERFORM INTEGRATION 383 

530 TN3»TN+DELTAT/3. 384 

CALL TBLU(NT1,0RD,TN3,TABQ,MQ) 385 

DYNP-ORD*QDELTA 386 

CALL TBLU(NT3,0RD,TN3,TABT,MT) 387 

THRUST=ORD»TDELTA 388 

CALL TBLU (NT3,ETM,TN3,ETMEAN,MEM ) 389 

CALL TBLU ( NT4, ETS, TN3, ETSIG, MES ) 390 

ET-ETM+RXDET*ETS 391 

AL- AUPR+THETN+AUD3l:TN3+THEDN:4:DELTAT/a . +THEDD*0 . 5*DELTAT**3 393 

TDDN3 = THEDDC+RATI01*ET*THRUST+4,775^:DYNP:tc(CNAS:»:AL+A3*AL*AL»AL) 393 
1*(CG-XCP0-DCPDA*ABS(AL) )/RINT 394 

THEDD“TDDN3 395 

THEDN1=>TDDN3:»:DELTAT+THEDN 396 

THETNl - (THEDNl+THEDN )2t:DELTAT/3.+THETN 397 

THEROl “THETNl-GDOT* ( TN+DELTAT ) 398 

XINNl»XINN+THEROR*DELTAT 399 

EPIT1-GINT*XINN1 400 

EPIN1-SIGN(AMIN1(ABS(EPIT1),XLIM),EPIT1 ) 401 

C TEST FOR START TIME FOR INTEGRAL GAIN CONTROL TERM 403 

IF(TN.LE.TKI)EPIN1»0. 403 

C TEST FOR DIUERGENT OUERSHOOT, THETA MAX,, ERROR MODE 3 404 

IF (ABS(THEROl)-TLIMIT ,GE, 0 ) GO TO 440 405 

IF (THEDDC .EG. 0 ) GO TO 550 -406 

C CONTROL MOTOR ON 407 

THETS-THER01+GAIN*THEDN1+EPIN1 408 



IF (ABS(THETS)-D .LT. 0 ) GO TO 560 

OUT OF DEADBAND. TEST FOR ONE INTEGRATION STEP UITHIN INTEGRAL 
GAIN START TINE 

IF ( ABS(TN-TKI)-DELT0 .GT. 0 ) GO TO 430 

UITHIN ONE STEP OF TKI. CHECK FOR INTEGRAL GAIN UALUE 

IF (ABS(EPINl) .EG. 0 ) GO TO 430 

IF (THETS .LT. O ) GO TO 540 

FC-FB 

THEDDC=TDDCB 
KCASE-0 
GO TO 600 
TOP NOTOR FIRING 
540 FC=FT 

THEDDC=TDDCT 
KCASE*0 
GO TO 600 

CONFUTE CONTROL ERROR SIGNAL 
CONTROL NOTOR OFF BUT CONING ON 
550 THETS=THER0H-GAIN*THEDN1+EPIN1 
TEST FOR OUERSHOOT 
560 IF (THEDN-GDOT) 570,590,580 
570 IF (THEDNl-GDOT .GE. 0 ) GO TO 590 
GO TO 600 

580 IF (THEDNl-GDOT .GT. 0 ) GO TO 600 
590 KCASE=1 

GO TO 250 
INCRENENT TINE 
600 TN=TN+DELTAT 

TINPN»TINPN+FC:*cDELTAT 

THEDN»THEDN1 

THETN=THETN1 

THEROR-THEROl 

XINN-XINNl 

TEST PRINT OPTION 




IF (NPRINT-a .LT. C ) GO TO 640 
IF (NPK .NE. 0 ) GO TO 6E0 
URITE( 6,870) 

URITEC 6,880) 

NLINE-NLINE+4 

NPK=NPK+1 

IF (NLINE+1-59 ,LE, 0 ) GO TO 630 
NPAGE-NPAGE+1 
CALL PAGEHD 
GO TO 610 

PRINT ATTITUDE ERROR, RATE, TINE AND IMPULSE 
URITE( 6,890) THEROR,THEDN,TN,TIMPN 
NLINE-NLINE+1 

IF (THEDDC ,EQ. 0 ) GO TO 650 
FC=0.0 
THEDDC=0,0 
D = DB 

GO TO 160 

SET DEADBAND FOR SUITCHING OFF 
D=DBOFF 

IF (THETS .LT. 0 ) GO TO 660 
FC*FB 

THEDDC*TDDCB 
GO TO 160 
TOP MOTOR FIRING 
FC-FT 

THEDDC»TDDCT 
GO TO 160 

RECORD FINAL UALUES FOR THIS CASE AND UPDATE SUMS 

THES( J )-THECAP 

TTCAP( J )-TCAP 

TIMPC J )-TIMPCP 

STS-STS+THECAP 

SSTS-SSTS+THECAP*THECAP 


43 

44 

45 

46 

47 

48 

49 

50 

51 
53 

53 

54 

55 

56 

57 

58 

59 

60 
61 
63 

63 

64 

65 

66 

67 

68 

69 

70 

71 
73 

73 

74 

75 

76 
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STT-STT+TCAP 477 

SSTT«SSTT+TCAP*TCAP 478 

STI-STI+TIMPCP 479 

SSTI-SSTI + TiriPCP*TinPCP 480 

IF (J-NE .LT. 0 ) GO TO 50 481 

COnPUTE THE MEAN AND STANDARD DEUIATION FOR SPECIFIED NUNBERdNC) 48E 
OF CAPTURE ANGLES, TIMES AND IMPULSES 483 

AM-NS 484 

BM=AM*(AM-1 .0) 485 

XMTS=STS/AM 48G 

XMTT=STT/AM 487 

XMTI*STI/AM 488 

SDTS=SQRT( (AM»SSTS-STS*STS)/BM) 489 

SDTT=SQRT( ( AM*SSTT-STT*STT )/BM ) 490 

SDTI=SQRT( < AM*SSTI-STI*STI )/BM ) 491 

C PRINT THE MEAN AND STANDARD DEUIATION AT SPECIFIED INTERUALS 49E 

IF (NPRINT .NE. 0 ) GO TO 690 493 

IF (NPH .GT. 0 ) GO TO 700 494 

690 NPAGE-NPAGE+1 495 

CALL PAGEHD 496 

URITE( 6,900) 497 

URITE( 6,910) 498 

NLINE«NLINE+4 499 

NPH'NPH+1 500 

700 IF (NLINE+1-59 .GT. 0 ) GO TO 690 501 

C PRINT RUNNING MEANS AND STANDARD DEUIATIONS 50E 

URITE< 6,9S0) N2,XMTS,SDTS,XMTT,SDTT,XMTI,SDTI 503 

NLINE-NLINE+1 504 

C CHECK FOR TOTAL NUMBER OF RUNS 505 

IF (J-NTRUN+INC . GE . 0 ) GO TO 710 506 

NE-NE+INC 507 

GO TO 50 508 

710 IF (J-NTRUN .GE. 0 ) GO TO 720 509 

NE-NTRUN 510 
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GO TO 50 511 

730 1-0 512 

C PRINT HISTOGRAM RESULTS 513 

IF (NHIST .EG, 0 ) GO TO 750 514 

C FILL HISTOGRAM CELLS UITH DATA 515 

DO 730 L-1,N2 516 

ICELL-(THES(L)-CMINH)/DIFC+1 517 

IF(ICELL.GT.NCEL)ICELL=NCEL 518 

CAHIST(ICELL)-CAHIST(ICELL)+1 519 

ICELL-(TTCAP(L)-TMINH)/DIFT+1 520 

IF(ICELL,GT.NCEL)ICELL-NCEL 521 

THIST(ICELL)=THIST(ICELL)+1 522 

ICELL-(TIMP(L)-PMINH)/'DIFP+1 523 

IF(ICELL.GT.NCEL)ICELL-NCEL 524 

730 HISTIP(ICELL)-HISTIP(ICELL)+1 525 

NPAGE=NPAGE+1 526 

CALL PAGEHD 527 


URITE< 6,1140) 528 

DO 740 L-1,NCEL 529 

CMINHl-CMINH 530 

TMINH1=TMINH 531 

PMINHl-PMINH 532 

CMINH=CMINH+DIFC 533 

TMINH=TMINH+DIFT 534 

PMINH-PMINH+DIFP 535 

PRINT OUT THE HISTOGRAM CELL DATA 536 

URITE( 6,1150 ) CniNHl,CMINH,CAHIST(L),TMINHl,TMINH,THIST(L),PMINH1537 
1,PMINH,HISTIP(L) 538 

740 CONTINUE 539 

GO TO 770 540 

750 NPAGE-NPAGE+1 541 

CALL PAGEHD 542 

URITEC 6,930) 543 

URITEC 6,940) 544 
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NLINE-NLINE+4 545 

760 I-I+l 546 

PRINT THE UNARRANGED ARRAVS OF CAPTURE ANGLE, TIME AND IMPULSE 547 

ONLY IF HISTOGRAM DATA IS NOT OUTPUT 548 

URITE( 6,950) THES ( I ) , TTCAP( I) ,TIf1P< I) 549 

NLINE=NLINE+1 550 

IF (N2-I ,LE. 0 ) GO TO 770 551 

IF (NLINE+1-59 .LE, 0 ) GO TO 760 55E 

GO TO 750 553 

770 CONTINUE 554 

C REARRANGE TOTAL IMPULSE, CAPTURE TIME AND ANGLE IN ASCENDING 0RDER555 

CALL ASCEND(N2,TIMP,0) 556 

CALL ASCEND(NE,TTCAP,0 ) 557 

CALL ASCEND(N2,THES,0) 558 

C PRINT THE REARRANGED ARRAYS IN ASCENDING ORDER 559 

1=0 560 

840 NPAGE=NPAGE+1 561 

CALL PAGEHD 562 

URITE( 6,960) 563 

URITEC 6,970) 564 

NLINE=NLINE+4 565 

850 1=1+1 566 

C COMPUTE DISCRETE PROBABILITY CC AND OUTPUT 567 

AA=N2+1 568 

BB-I 569 

CC=BB/AA 570 

URITE( 6,980) CC , THES ( I ) , TTCAP ( I ) , TIMP ( I ) 571 

NLINE-NLINE+1 572 

IF (N2-I .LE. 0 ) GO TO 860 573 

IF (NLINE+1-59 .LE. 0 ) GO TO 850 574 

GO TO 840 575 

860 CONTINUE 576 

C PRINT THE TOTAL NUMBER OF DIUERGENT CASES ENCOUNTERED 577 

NPAGE-NPAGE+1 578 
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CALL PAGEHD 
URITE( 6,990) NERR3 
URITE( 6,1000) NERR1,NERR4 
URITE( 6,1010) K 
GO TO 10 


870 rORriAT (//59H 
lULSE ) 

880 FORflAT (59H 
IC)/ ) 

890 FORMAT 
900 FORMAT 
IMA 

910 FORMAT 


THETA 


(DEG) 


THETADOT 


(DEG/SEC) 


TIME 


(SEC) 


(5X,F10.4,3(F14.4 ) ) 
(//79H ITER 

MEAN SIGMA 

(79H NO 


MEAN 


SIGMA 


MEAN 


) 


THETA 


THETA 


IMPULSE 


IMPULSE /) 


920 FORMAT 
930 FORMAT 
940 FORMAT 
950 FORMAT 
960 FORMAT 
ILSE ) 
970 FORMAT 
1 ) /) 

980 FORMAT 
990 FORMAT 
1 

1000 FORMAT 
1 
2 


(5X,I6,F13.5,5(F11.5) ) 
(//47H CAP. ANGLE 

(47H (DEGREES) 

(5X,F11.5,F15.5,F16.5) 
(//SSH PROB 


CAP. TIME 
(SECONDS) 


TIME 


IMPULSE 

(LB-SEC)/) 


(57H 


CAP. ANGLE 


(DEGREES) 


CAP. TIME 


(SECONDS) 


DIUERGENT CAPTURES 


(Fll .3,F15.5,F15.5,F16.5) 

(//, 1X,37HT0TAL NUMBER OF 
15) 

(///,1X,41HT0TAL NUMBER OF ABNORMAL CASES 
42HDEADBAND OUERSHOOT 

42HTIME LIMIT UITHOUT DIUERGENCE OR CAPTURE - 


DELETED: 

,I5,/,5X, 

^ 15 ,/) 


1010 FORMAT (//52H 
15X, 115 ) 

1020 FORMAT ( 14(15 ) ) 
1030 FORMAT (13A6) 


FINAL UALUE OF RANDOM SEQUENCE INTEGER 


579 

580 

581 

582 

583 

584 
IMP585 

586 

(LB-SE587 

588 

589 
SIG590 

591 
TIME592 

593 

594 

) 595 

596 

597 
IMPU598 

599 

(LB-SEC600 

601 

602 

603 

604 

, // , 5X, 605 
606 
607 
,608 

609 

610 
611 
612 
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1040 FORMAT (110) 613 

1050 FORMAT (aCElS.S)) 614 

1060 FORMAT (E15.8) 615 

1070 FORMAT ( I5/, ( 6E10 . 3 ) ) 616 

1100 FORMAT (/,10X,10HERROR MODE , IE, 14HENC0UNTERED - ,3A6,1X, 617 

1 aOHUITH FOLLOUING DATAt , lEX, 1 IHTIME (SEC )-,F7. 1,3X, 618 

a 14HRATE(DEG/SEO- , F7 . E, 3X, 16HATT . ERROR( DEG ) - , F8 ♦ 3, 619 

3 iaX,llHO (PSF) » ,F7.1,3X,14HALPHA(DEGREE)» >F7.a,3X, 6B0 

4 16HTHRUST MIS ( DEG ) = ^ F8 . 3, /, lEX, IIHF TOP(LB) - ,F7,1,3X, 6B1 

5 14HF LOUER (LB) - , F7 . B, 3X, 16HLAST OOERSHOOT -,F8.3,/) 6BB 

1140 FORMAT (1H1,4X,B0HCAPTURE ANGLE NO ,8X,4HTIME,7X,BHN0,7X, 6B3 

1 7HIMPULSE,10X,BHNO,/) 6B4 

1150 FORMAT (1X,F7.B,3H TO>F10.B, IH, , 14, F7.B,3H T0,F6.B, IH, , I4,F9.B, 6B5 

1 3H T0,F9,a,lH,,I4) 6B6 

END 6B7 
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> 
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*DECK ASCEND 628 

SUBROUTINE ASCEND (L,UAL,I1) 629 

SUBROUTINE FOR REARRANGING ARRAV IN ASCENDING ORDER 630 

OR DESCENDING ORDER IF N-1 631 

DIMENSION UAL(l) 632 

K-L-1 633 

DO 40 634 

KB-J+1 635 

DO 40 JL=»KB,L 636 

IF(M) 10,10,20 637 

10 IF (UAL( J )-UAL( JL) ,LE. 0 ) GO TO 40 638 

GO TO 30 639 

20 IF(UAL( JL)-UAL( J ).LE.0) GO TO 40 640 

30 TENP-UAL(J) 641 

UAL(J)-UAL(JL) 642 

UAL(JL)*TEMP 643 

40 CONTINUE 644 

RETURN 645 

END 646 



*DECK PAGEHD 647 

SUBROUTINE PAGEHD 648 

C THIS SUBROUTINE PRINTS HEADING AT THE TOP OF EACH PAGE 649 

COMMON RXD,Q,K,NPAGE,NLINE,NRUN,NTITLE(a6) 650 

NLINE-8 651 

URITE( 6^10) 65a 

URITE( 6>a0) NRUN^NPAGE 653 

URITE(6,5) (NTITLE(J)>J-l>a6) 654 

RETURN 655 

5 FORMATC/, ( IX, 13A6 ) ) 656 

10 FORMAT (IHl) 657 

E0 FORMAT <5X,7HRUN NO . , 15, 49X, 8HPAGE NO., 15) 658 

END 659 
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*DECK RNDX 660 

SUBROUTINE RNDX(L) 661 

THIS SUBROUTINE COMPUTES THE PSEUDO-RANDOM NUMBER 663 

L-0 GIUES RANDOM NORMAL DEUIATE RXD FROM N(0>1) 663 

L-1 GIUES RANDOM UNIFORM NUMBER Q 0 TO 1 U(0,1) 664 

COMMON RXD,Q,K,NPAGE,NLINE,NRUN,NTITLE(a6) 665 

DATA A/1«0E07/, JJ/10000000/ 666 

KA-KXE3 667 

N-KA/JJ 668 

M-N^:JJ 669 

K=KA-M-N 670 

AK=K 671 

C COMPUTE RANDOM NUMBER Q UNIFORM FRON ZERO TO ONE 672 

Q=AK/A 673 

IF ( L . GT . 0 ) RETURN 674 

C CHANGE Q TO PLUS AND MINUS 0 TO 0.5 RANGE 675 

IF (Q-.5 .LE. 0 ) GO TO 10 676 

> Q=l.-Q 677 

AB--1.0 678 

^ GO TO 20 679 

C 680 

10 AB=1.0 681 

C TRANSFORM Q TO GAUSSIAN RANDOM NORMAL DEUIATE RXD 682 

20 Y=SQRT(-ALOG(Q*Q ) ) 683 

D=2.515517+0.802a53XY+0.010328XV*Y 684 

E = 1 .0+1 .432788*Y+0. 1892692tcV*Y+0.001308*YXY*V 685 

RXD*ABX(Y-(D/E ) ) 686 

RETURN 687 

END 688 

689 
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SUBROUTINE TBLU (NT,Y,X,T,P1) 690 

C SINGLE TABLE LOOKUP SUBROUTINE 691 

C NT - NUMBER OF UALUES IN ARRAY 692 

C Y - RETURNED ORDINATE 693 

C X - ABSCISSA UALUE CALLED 694 

C T “ INPUT TABLE OF ALTERNATING ABSCISSAS AND ORDINATES 695 

C ORDINATES MUST BE MOHOTONICALLY INCREASING 696 

C M = PREUIOUS INDEX USED IN THIS TABLE LOOKUP 697 

C THIS INDEX GETS CHAN .ID TO CURRENT UALUE 698 

DIMENSION T(1 ) 699 

10 IF (T(M)-X) 50,20,30 ^ 700 

20 Y = T(M+1 ) 701 

RETURN 702 

30 IF (T(l)-X ,LT. 0 ) GO TO iO 703 

M-1 704 

GO TO 20 705 

40 M=n-2 706 

GO TO 10 707 

50 MM=M+2 708 

IF (MM-NT-1 .LE. 0 ) GO TO 60 709 

M“NT-1 710 

GO TO 20 711 

60 IF (T(MM)-X .GT. 0 ) GO TO 70 712 

M=MM 713 

GO TO 50 714 

70 M=MM-2 715 

DT = T(MM )-T(M ) 716 

IF (DT .NE. 0 ) GO TO 80 717 

Y=T(M+1) 718 

RETURN 719 

80 DY = T(MM+1 )-T(M+l ) 720 

DDT-X-T(M) 721 

Y-T(M + 1 ) + DY:*:DDT/'DT 722 

RETURN 723 

END 
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